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ABSTRACT 

This paper presents an analysis of the correlation between the number of globular clusters (iVcc) 
in giant galaxies and the mass of the galaxies' central supermassive black hole (AfsMBii)- I construct 
a sample of 20 elliptical, spiral, and SO galaxies with known SMBH masses and with accurately- 
measured globular cluster system properties derived from wide-field imaging studies. The coefficients 
of the best- fitting Nqc— -Msmbh relation for the early- type galaxies are consistent with those from 
previous work but in some cases have smaller relative errors. I examine the correlation between Aqc 
and Msmbh for various subsamples and find that elliptical galaxies show the strongest correlation 
while SO and pseudobulge galaxies exhibit increased scatter. I also compare the quality of the fit 
of the numbers of metal-poor globular clusters versus SMBH mass and the corresponding fit for 
metal-rich globular clusters. I supplement the 20-galaxy sample with ten additional galaxies with 
reliable iVcc determinations but without measured AfsMBH- I use this larger sample to investigate 
correlations between Nqq and host galaxy properties like total galaxy luminosity and stellar mass and 
bulge luminosity and mass. I find that the tightest correlation is between Nqc and total galaxy stellar 
mass. This lends support to the notion that Aqc and Msmbh are not directly linked but are correlated 
because both quantities depend on the host galaxy potential. Finally, I use the Nqc— -Msmbh relation 
derived from the 20-galaxy sample to calculate predicted Msmbh values for the ten galaxies with 
accurate Agc measurements but without measured SMBH masses. 

Subject headings: black hole physics - galaxies: elliptical and lenticular, cD - galaxies: spiral - galaxies: 
star clusters - galaxies: formation - globular clusters: general 



1. INTRODUCTION 

It is now well-known that a correlation exists between 
the mass of the supermassive black hole (SMBH) at the 
centers of giant galaxies and the velocity dispersion of 
the galaxies' spheroidal component (Ferrarese & Mer- 
ritt 2000, Gebhardt et al. 2000b and many subsequent 
papers). This u M-a relation" has been interpreted as 
indicating a close connection between the formation and 
evolution of galaxy bulges and central SMBHs. Much 
observational and theoretical work has been devoted to 
establishing the exact form of the relation (e.g., Tremaine 
et al. 2002, Gultekin et al. 2009) as well as the reasons 
for its existence (e.g., Burkert & Silk 2001, Miralda- 
Escude & Kollmeier 2005, Di Matteo et al. 2005, Robert- 
son et al. 2006, Johansson et al. 2009, Jahnke & Mac- 
cio 2011). The general consensus seems to be that the 
presence of SMBHs at the centers of giant galaxies is 
a natural consequence of hierarchical structure forma- 
tion and fits well within our overall picture of galaxy 
assembly (e.g., Robertson et al. 2006, Jahnke & Maccio 
2011). We also know that SMBHs must have formed rel- 
atively early in many galaxies' histories: observations 
of high-redshift quasars demonstrate that SMBHs ex- 
isted at z of ^6—7 (e.g., Fan et al. 2006, Mortlock et 
al. 2011), when the Universe was less than a Gigayear 
old. On the other hand, the detailed physical mecha- 
nisms responsible for seeding the central SMBHs in the 
first place, and then growing them over time, are not yet 
well- understood (e.g., Volonteri & Rees 2006, Omukai, 
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Schneider, & Haiman 2008, Mayer et al. 2010). 

Like SMBHs, many globular clusters (GCs) appar- 
ently also formed during the early stages of galaxy for- 
mation. GCs are luminous (My —11 to —4), compact 
(i?i/ 2 of a few parsecs), populous (~10 4 — 10 6 members) 
star clusters that orbit their host galaxies at galacto- 
centric distances of less than a kiloparsec out to hun- 
dreds of kiloparsecs (e.g., Rhode & Zepf 2004, Rhode et 
al. 2007). GCs in the Milky Way, which arguably have 
the most accurate absolute age determinations, typically 
have ages of ~11 — 13 Gyr (e.g., Dotter et al. 2010, Forbes 
& Bridges 2010, and references therein). Observations of 
external galaxies have detected GCs with intermediate 
ages (e.g., 1.5—4 Gyr; Goudfrooij et al. 2007) as well as 
young (<1 Gyr) star clusters, with masses equal to or 
greater than GCs (e.g., Bastian et al. 2006, Whitmorc 
et al. 2010), that were likely produced in gas-rich galaxy 
mergers. Because they appear to have formed during in- 
tense star formation events, including events triggered by 
mergers, the properties of GCs can be used to trace the 
major assembly and evolutionary episodes of their host 
galaxies (Ashman & Zepf 1998; Brodie & Strader 2006). 
The total numbers, spatial distributions, ages, metallic- 
ities, and kinematics of GCs provide important physical 
clues regarding the origin and evolution of their parent 
galaxies. 

A few recent papers have begun to explore possible 
links between these two constituents of galaxies, SMBHs 
and GCs. Spitler & Forbes (2009) presented a method 
for using the total number of GCs (Aqc) in a galaxy 
to estimate the mass of the galaxy's halo. They in turn 
showed that, for a sample of about a dozen galaxies, the 
halo mass derived from the number of GCs correlates 
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fairly well with the measured SMBH mass. Burkert & 
Tremaine (2010; hereafter BT10) investigated a more di- 
rect correlation between SMBHs and GCs by showing 
that there is a tight correlation between Nqq and the 
mass of the central SMBH (Msmbh) in early- type galax- 
ies. Using a sample of 13 galaxies with Nqc estimates 
and measurements of SMBH mass, BT10 found that the 
relationship between Nqq and Msmbh was even tighter 
than the M-a relation for those same galaxies. BT10 
also found that the mass of the central SMBH in the el- 
liptical and SO galaxies they studied was roughly equal 
to the mass in GCs. They concluded that the origin 
of the Nqq — Msmbh correlation was obscure and that a 
larger sample of galaxies was needed to further explore 
its causes and implications. 

Harris & Harris (2011; hereafter HH11) followed up on 
this previous work by expanding the sample of galaxies 
to more than twice that of BT10. Starting with the list 
of galaxies for which Msmbh had been measured, they 
selected from the literature those with either published 
Nqq values or with information about the GC system 
that would allow them to produce at least rough esti- 
mates of Nqq. (BT10 had restricted their sample to 
galaxies included in the ACS Virgo Cluster Survey pa- 
per on GC populations by Peng et al. 2008 or the list of 
GC system properties compiled by Spitler et al. 2008.) 
The end result was a sample of 33 elliptical, SO, and spi- 
ral galaxies with both Msmbh and Nqq values. HH11 
used the elliptical galaxies in this sample to derive a re- 
lation between Nqc and Msmbh and found that it was 
consistent with the relation from BT10 within the errors. 
They also examined how the correlation varied with the 
morphology of the host galaxies. HH11 reasoned that 
correlations between GC populations and SMBHs may 
arise because of similarities in the conditions and epoch 
of formation of these two types of objects. They noted 
that both likely originate in situations in which high gas 
densities are produced by energetic collisions and merg- 
ers of gas clouds. Such circumstances would be expected 
to occur in the early universe (e.g., z ~10— 15 or higher), 
during the initial phase of galaxy assembly and forma- 
tion. 

Snyder, Hopkins, & Hcrnquist (2011) drew on the data 
sets and results of both BT10 and HH11 in order to fur- 
ther examine the possible causes of the Nqq— Msmbh 
correlation. They argued that the observed correlation 
is a result of a link between SMBH mass and the binding 
energy of the host galaxy bulge - the so-called Black Hole 
Fundamental Plane (BHFP; e.g., Hopkins et al. 2007) - 
combined with a link between Nqq and the stellar mass 
of galaxy bulges at fixed velocity dispersion. They ex- 
amined in detail the residuals in the Nqq — Msmbh cor- 
relation as well as correlations between Nqq and other 
galaxy properties. Snyder et al. showed that although 
the scatter in the Nqq — Msmbh correlation is small, this 
does not imply a "special" relationship betwen GCs and 
SMBHs, but instead is due to both Msmbh and Nqq 
correlating with bulge mass. They point out that while 
the relationship betwen SMBH mass and bulge potential 
is consistent with our general understanding of galaxy 
mergers and physical processes like black hole accretion, 
feedback, and gas cooling, the reasons why Nqq should 
be so tightly correlated with bulge stellar mass are not 
entirely clear. 



The measurements of Nqq used in these papers are 
taken from heterogeneous data sets of varying quality. 
Therefore a key question that should be explored is what 
would happen to the relation between Nqc and Msmbh, 
as well as the relations between Nqc and other fun- 
damental galaxy properties, if the galaxy sample were 
constrained to include only galaxies with the most well- 
determined GC system properties. Indeed, BT10 con- 
clude their paper by noting that, "An important next 
step is to expand the sample of galaxies having both re- 
liable SMBH masses and reliable GC populations." I am 
leading an ongoing wide-field CCD imaging survey aimed 
at establishing the total numbers and other global prop- 
erties of the globular cluster systems of giant galaxies 
(Rhode & Zepf 2001, 2003, 2004; Rhode et al. 2005, 2007, 
2010; Hargis et al. 2011; Young, Dowell, & Rhode 2012, 
Hargis & Rhode 2012). Using the survey results to inves- 
tigate the issues raised by BT10, HH11, and Snyder et al. 
(2011) is a natural application of the data. Accordingly, I 
have compiled a sample of 20 galaxies with measured GC 
system properties from observational studies that meet 
specific criteria - i.e., studies in which the imaging data 
cover a fair fraction of the GC system and in which quan- 
tities like the global color fraction are known. To date, 
we have observed ~25 galaxies for the wide- field imaging 
survey; eight of these have SMBH mass determinations 
in the literature, so they make up eight of the 20 galax- 
ies included in this paper. The remaining twelve galaxies 
included here are drawn from several GC system studies 
in the literature. 

Here I present the results of this investigation. In Sec- 
tion 2, I discuss how Nqq and other GC system proper- 
ties are measured in giant galaxies, what issues give rise 
to uncertainties in those numbers, and the criteria I ap- 
plied to decide what data to include in the current study. 
I then present the data set, including details about how 
measurements were derived for each individual galaxy. 
Section 3 describes the results of the study: I examine 
correlations between Nqq and SMBH mass for the sam- 
ple galaxies, investigate the link between metal-poor and 
metal-rich GCs and SMBH mass, and look for trends 
in the Nqq— Msmbh correlation with galaxy morphol- 
ogy and bulge type. I also investigate the link between 
the numbers of GCs and host galaxy properties. In the 
last subsection of Section 3, I calculate predicted SMBH 
masses for a set of galaxies from my wide-field survey. 
The final section of the paper gives a brief summary of 
the main conclusions. 

2. THE DATA SET 

2.1. Deriving Global Properties of GC Systems 

Quantifying the properties of the GC systems of gi- 
ant galaxies beyond the Local Group is a challenging 
task. Although it has been known for decades that gi- 
ant galaxies besides the Milky Way and M31 host GCs 
- e.g., authors like Baum (1955), Sandage (1961), and 
Dawe & Dickens (1976) speculated that the point-like 
objects surrounding galaxies in the Virgo and Fornax 
Clusters were GCs - it wasn't until the early 1980s that 
systematic studies of GC populations in large numbers 
of galaxies began in earnest (Harris 1991). GCs in galax- 
ies more than a few Mpc away are unresolved in typical 
(~0.5— 1.0") ground-based seeing; the median half-light 
radius of a Milky Way globular cluster is 3 pc, which 
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translates to —0.04" at the distance of the Virgo Cluster 
(Ashman & Zepf 1998). The luminosities of GCs follow a 
roughly Gaussian distribution referred to as the GC Lu- 
minosity Function (GCLF), which has a peak magnitude 
M£ = -7.4±0.25 and a dispersion a = 1.4±0.2 mag (e.g., 
Whitmore et al. 1995, Kavelaars et al. 2000, Barmby et 
al. 2001, Kundu & Whitmore 2001, Jordan et al. 2007). 
Therefore, studies of extragalactic GC systems aim to 
identify GCs as a population of point sources with optical 
magnitudes and colors like GCs arrayed in a centrally- 
concentrated distribution around the host galaxy. To 
derive the total number of GCs in a galaxy's system, one 
first identifies and counts the GCs (trying to distinguish 
them from other compact objects like foreground stars or 
background galaxies) and then makes various corrections 
to account for missing spatial coverage, the detection lim- 
its of the images, and contamination from non-GCs, to 
arrive at a final estimate of Nqc- 

Early studies of extragalactic GC systems employed 
photographic plates and source detection in one or two 
filters to identify GCs and separate them from contam- 
inating objects. GCs can be found tens to hundreds of 
kiloparsecs from the center of their host galaxy (Har- 
ris 1991; Ashman & Zepf 1998) and large-format photo- 
graphic plates were a good match for the extended na- 
ture of giant galaxy GC systems. On the other hand, 
the low quantum efficiency of plates meant that only 
the brightest portion of the GC population could be de- 
tected. Photographic studies also suffered from high lev- 
els of contamination from stars and galaxies. The result 
was that quantities like Ngc and color distributions for 
the full system were often highly uncertain and based 
on observations of only a small percentage of the total 
GC population. A few examples arbitrarily selected from 
the review article by Harris (1991) serve to illustrate this. 
Harris summarized the progress of extragalactic GC sys- 
tem studies and compiled estimates of Nqc from the lit- 
erature, listing both number of clusters observed (N f, s ) 
and the number of GCs in the system. The latter quan- 
tity is computed by extrapolating N Q b s over all magni- 
tudes and radii. For the Sab galaxy NGC 4569, N b s 
= 30±10 and the derived number of GCs is 1000±400. 
For the SO galaxy NGC 3607, N obs = 50±35 and derived 
A GC = 800±600; for the elliptical NGC 3311, N obs = 
414±31 and the estimated Nqc is 41 times larger, at 
17000±6000. 

The increased availability of CCD detectors in the 
late 1980s made detecting faint GCs in external galax- 
ies more efficient and made it easier to use color cri- 
teria (e.g., B — V) to select GCs, thereby reducing con- 
tamination from Galactic stars and background galaxies. 
Early CCDs had small formats, however, which meant 
that large radial extrapolations were necessary to de- 
rive global GC system properties. For instance, a careful 
imaging study by Lee, Kim, & Geisler (1998), performed 
using a 2048 x 2048-pixel CCD and Washington C and 
Tl filters, traced the GC system of the Virgo elliptical 
NGC 4472 (M49) to 7' (-34 kpc). Lee et al. calculated 
both Nqc and specific frequency S N for NGC 4472's GC 
system. Specific frequency was introduced by Harris & 
van den Bergh (1981) and is the number of GCs nor- 
malized by the host galaxy V-band luminosity: 

S N =N GC 10 +0 - 4(Mv+15) (1) 



To calculate total number and Sn, Lee et al. integrated 
their observed GC system radial distribution to 10'. 
Later BVR imaging with an 8192 x 8192-pixel mosaic 
CCD camera showed that NGC 4472's GC system actu- 
ally extends to at least 23', or > 110 kpc (Rhode & Zepf 
2001). 

Beginning in the early 1990s, the superior resolution 
and sensitivity of the Hubble Space Telescope (HST) 
cameras enabled much deeper imaging of extragalactic 
GC systems as well as allowing faint background galaxies 
to be resolved and therefore eliminated from GC studies. 
In many ways these advances revolutionized the study of 
GC populations beyond the Local Group, allowing the 
GCLF, GC color distributions, and even GC sizes to be 
studied in detail in many galaxies (see Brodie & Stradcr 
2006 and references therein). Even so, HST's small ficld- 
of-view (FOV) meant that typically only a fraction of a 
galaxy's GC system was imaged, making it difficult or 
impossible to quantify global quantities like total num- 
ber and color distribution for the system. For example, 
the HST ACS Virgo Cluster Survey (ACSVCS; Cote et 
al. 2004) imaged 100 Virgo ellipticals and SO galaxies 
to radii of 2.4' (12 kpc), whereas early-type galaxy GC 
systems at Virgo Cluster distances often extend an order 
of magnitude beyond that, to 10'— 20', or —50—100 kpc 
(Rhode & Zepf 2001, 2004). To arrive at estimates of 
Aqc and global Sn for galaxies from the ACSVCS, Peng 
et al. (2008) supplemented their survey data with paral- 
lel WFPC2 imaging (that provided observations of the 
outer regions of the GC systems of massive galaxies) 
and ground-based surface density profiles from wide-field 
imaging studies in the literature. Some of the Agc and 
Sn values derived from the ACSVCS end up being quite 
similar to previously-published values, while others are 
significantly different. For example, the Harris (1991) lit- 
erature compilation lists the number of GCs in the Virgo 
E6 galaxy NGC 4564 as A GC = 1200±400; Peng et al. es- 
timate Aqc — 213±31 for this galaxy. Harris (1991) lists 
A GC = 2200±440 for the Virgo E5 galaxy NGC 4621 and 
Peng et al. derive A GC = 803±305. The Harris (1991) 
table lists A GC = 3500±1200 for the Virgo E2 galaxy 
NGC 4365 and Peng et al. measure a similar number, 
although with a much smaller error: Nqc — 3246±598. 
The variation in these published values serves both to 
highlight the very real difficulty of measuring global GC 
system properties (especially Nqq and Sn) and to mo- 
tivate why it is important to critically examine how the 
numbers were derived before using them to study other 
properties of the parent galaxies. 

In the late 1990s, large-format and mosaic CCD im- 
agers became available on many 4-meter-class telescopes, 
making it possible to image much larger portions of gi- 
ant galaxy GC systems outside the Local Group in one 
or a few pointings. Because of the improved efficiency of 
this approach — both in terms of increased sensitivity 
compared to photographic studies, and increased areal 
coverage compared to small-format CCDs or HST imag- 
ing — it also became more common to image the GC 
populations in more than two broadband filters (e.g., 
Rhode & Zepf 2001, Tamura et al. 2006, Spitlcr et al. 
2008). The observational goals of these types of studies 
typically were to cover the majority of the spatial extent 
of the GC systems, to do deep photometry in order to 
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sample a significant fraction of the GC luminosity func- 
tion (GCLF), and to reduce contamination from stars 
and galaxies by selecting GCs via their magnitudes and 
colors in multiple filters. The scientific objectives were to 
accurately measure global values for Agc, Sn, and the 
distributions of GC colors (which should vary primarily 
due to the metallicities of the clusters, with metal-poor 
GCs having blue broadband colors and metal-rich GCs 
appearing red) as well as to examine how these quantities 
change with galactocentric radius. As Brodie & Strader 
(2006) note in their review article on extragalactic GC 
systems, relatively few of these types of studies have been 
done to date, but the pattern so far is that the Nqc and 
Sn values they produce are often smaller than past val- 
ues and the errors on A^gc an d Sn are typically reduced 
by a factor of two or more (e.g., Rhode & Zcpf 2001, 
2003, 2004, Spitler et al. 2008). For example, Gomez & 
Richtler (2004) used the Calar Alto 3.5-m telescope and 
MOSCA CCD detector to image three pointings around 
the Virgo galaxy NGC 4374. They observed the GC sys- 
tem radial profile out to 12' (~60 kpc) and derived Nqc 
= 1775±150 and Sn = 1.6±0.3, whereas the previous 
value was N GC = 3400±800 and S N = 5.6±1.3 from a 
photographic study by Hanes (1977). 

It is likely that a number of factors contribute to the 
smaller total numbers and specific frequencies derived 
from modern wide-field CCD studies. One reason is that 
contamination from non-GCs is significantly reduced in 
modern studies because GCs are being selected according 
to their magnitudes in two or more filters and because 
good image resolution (<1") allows many background 
galaxies to be resolved and discarded from GC lists. An- 
other contributing factor is that power laws are often 
used in small-field imaging studies to fit the radial pro- 
file of a GC system. The power law, of the form log a G c 
= aO + al log r, where r is projected radius and ctqc 
is the surface density of GCs, is then integrated over all 
projected radii, or out to some assumed value for the ra- 
dial extent of the system. We have found in our ongoing 
GC system survey that de Vaucouleurs law profiles of the 
form log (tqq = aO + al r 1 / 4 arc often a better fit for 
GC system profiles out to large radius, and these drop 
off more quickly in the outer regions than a power law 
profile would (e.g., Rhode & Zepf 2004). Furthermore, 
we integrate the profiles only out to the radius where 
the GC system surface density falls off to zero within the 
errors, rather than integrating over all radii or to some 
arbitrary value. Both of these issues likely contribute to 
reduced Agc an d Sn values and uncertainties derived 
from wide-field studies that trace the full extent of the 
GC system. In some cases the newer studies have deeper 
imaging that enables a better determination of what frac- 
tion of the GCLF has been imaged, which in turn yields a 
more accurate correction for magnitude incompleteness. 
Finally, in a few cases - for example, a 2003 study of the 
GC system of the Fornax galaxy NGC 1399 by Dirsch et 
al. — a smaller or more accurate GC specific frequency 
may be derived simply because the total galaxy mag- 
nitude was revised based on improved wide-field CCD 
surface photometry. 

Finally, it is worthwhile to briefly note here the effect 
that errors in measured galaxy distances have on deter- 
minations of Agc and S n for extragalactic GC systems. 



Authors of GC system studies typically adopt a distance 
to the host galaxy from the literature. This distance is 
then folded into the calculations of distance-dependent 
quantities like the total galaxy magnitude and the frac- 
tional coverage of the theoretical GCLF. The latter quan- 
tity is usually determined by (1) assuming a value for the 
peak absolute magnitude and dispersion of a Gaussian 
GCLF and (2) fitting the observed luminosity function 
of GC candidates to the Gaussian by varying the normal- 
ization and then calculating how much of the area under 
this theoretical curve has been covered by the data. The 
absolute peak magnitude and/or dispersion of the theo- 
retical GCLF are sometimes also varied by a few tenths 
of a magnitude. Typically Agc and Sn determinations 
in the literature include errors on the counts of GC can- 
didates and contaminants and perhaps some modest un- 
certainty associated with the GCLF fitting process, but 
do not explicitly include the uncertainties in the galaxy 
distance (see, e.g., Harris 1991, Ashman & Zepf 1998, 
Spitler et al. 2008 and other compilations of GC sys- 
tem properties). If the distance to a galaxy has been 
underestimated and the true distance is greater than as- 
sumed, the calculated fractional GCLF coverage will be 
larger than it should be and the final Agc will be an 
underestimate of the true value. (That is, if the distance 
were corrected to its larger value, the final A^gc would 
increase.) Because both the galaxy magnitude and Agc 
contribute to Sn, distance errors have a different net ef- 
fect on specific frequency estimates. A galaxy that in re- 
ality is farther away than the adopted distance will have 
an underestimated Nqc > but its intrinsic luminosity will 
likewise be underestimated. Changes in Agc and My in 
Equation 1 will counteract each other, so an increase in 
distance to the galaxy can produce a larger Agc but a 
smaller GC specific frequency (see an example of this in 
our study of the GC system of NGC 7814; Rhode & Zepf 
2003). 

2.2. Constructing the Sample 
2.2.1. Criteria Used to Select Galaxies 

Given the issues described in the previous section, my 
objective was to compile a sample of giant galaxies with 
well-determined measurements of A^gc and GC specific 
frequency, drawing from my own survey and from the 
literature. I also decided to restrict the sample to those 
galaxies for which the global GC color distribution (i.e., 
the number of GCs versus broadband color, which is an 
indicator of metallicity) is known. Model scenarios for 
the formation and evolution of GC systems (e.g., Ash- 
man & Zepf 1992, Cote at el. 1998, Forbes et al. 1997, 
Beasley et al. 2002, Muratov & Gnedin 2010) typically 
make predictions for the ratio of blue (metal-poor) and 
red (metal-rich) GCs and how that ratio should change 
with radius in different types of galaxies. Because the 
metal-rich GC population has been shown to be more 
centrally-concentrated than the metal-poor population in 
some galaxies (e.g., Lee et al. 1998, Rhode & Zepf 2004), 
color distributions measured in only the central portions 
of GC systems may not accurately represent the global 
color distribution; thus wide-field coverage can also be 
important for measuring the true global fractions of blue 
and red GCs in a galaxy. 

Rhode, Zcpf, & Santos (2005) examined how the spe- 
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cific frequencies of blue, metal-poor GCs varied with host 
galaxy stellar mass for a sample of giant spiral, SO, and 
elliptical galaxies. We formulated a set of criteria that 
we believed would limit the sample to well-determined 
measurements. I have adopted the same standards here. 
The GC system studies included in the current sample 
must meet the following criteria: (1) at least 50% of the 
estimated radial extent of the GC system must have been 
observed; (2) imaging data must have been acquired in 
at least two filters, so that contamination from non-GCs 
can be reduced and the GC color distribution can be 
quantified; (3) an estimate of Nqq is given or can be de- 
rived in a straightforward way from the published data, 
and (4) the 1-er error on Nqc ot &n must be <40%. 
Since my aim for the current study is to use the galaxy 
sample to investigate the connection between supermas- 
sive black holes and GC populations, the galaxies in the 
sample must also have a measurement of SMBH mass 
in the literature. I used the lists of galaxies in BT10 
and HH11 for initial guidance regarding which galaxies 
had SMBH mass measurements, but also looked through 
the literature for additional measurements that may not 
have appeared in either of those papers. The papers 
from which I drew most of the SMBH mass measure- 
ments are the compilations of Gultckin et al. (2009) and 
Graham (2008) and references therein. I also examined 
the updated compilation of SMBH masses in Graham et 
al. (2011), but found no additional Msmbh values that 
would supplement the numbers taken from the Gultckin 
et al. (2009) and Graham (2008) papers. 

The final result of my search is a list of 20 galaxies 
that meet the stated criteria for GC system observations 
and also have published SMBH mass measurements. The 
data for these galaxies and the references from which the 
quantities are derived are given in Table 1. The table 
lists, in this order: galaxy name and morphological type; 
SMBH mass and associated reference(s); velocity disper- 
sion (for the galaxy or bulge, as appropriate); galaxy ab- 
solute V magnitude; the galaxy distance (in Mpc) that 
was assumed for the GC system values; number of GCs; 
GC specific frequencies; fraction of blue GCs in the sys- 
tem; and the reference for the GC system properties. 
(Note that estimated uncertainties on the fraction of blue 
or red GCs in the galaxies in this sample are typically 
a few to <10%; Rhode & Zepf 2001, 2004, Peng et al. 
2006). The GC specific frequency Sn in column 9 of Ta- 
ble 1 is as defined in Equation 1. Another type of specific 
frequency, T, is given in column 10. T was suggested by 
Zepf & Ashman (1993) and is defined as: 
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where Nqq is the number of GCs and Mq is the stellar 
mass of the host galaxy. Zepf & Ashman (1993) point 
out that Sn is affected by variations in V-band stellar 
mass-to-light ratios for galaxies of different morpholog- 
ical types and stellar populations, so T can be useful 
when one is comparing GC specific frequencies for galax- 
ies over a range of morphologies. I have adopted the 
mass-to-light ratios used by Zepf & Ashman when they 
defined T: M/L v = 10 for ellipticals, M/L v = 7.6 for SO 
galaxies, and M/Ly from 6.1 to 4.0 for spiral galaxies, 
with M/Ly decreasing with later Hubble type. 



2.2.2. Comments on Individual Galaxies 

Galaxies from Our Wide-Field GC System Survey. — 
The GC system values for eight of the galaxies in Ta- 
ble 1 come from the wide-field GC system survey that 
I have been leading (Rhode & Zepf 2001, 2003, 2004; 
Rhode et al. 2005, 2007, 2010; Hargis et al. 2011; Young 
et al. 2012, Hargis & Rhode 2012). We image spiral, SO, 
and elliptical galaxies at distances of ~10— 25 Mpc with 
large-format and mosaic CCD imagers. The FOVs of the 
cameras we use are ~10— 36' on a side, which translates 
to ^30—200 kpc at these distances; this is sufficient to 
observe the full radial extent of the GC systems of the 
target galaxies, which typically range from ~ 10—100 kpc 
(e.g., Rhode et al. 2010). We select the GC candidates 
around each galaxy via three-color photometry to both 
reduce contamination from non-GCs and to allow us to 
investigate the color distributions and color gradients of 
the GC systems. The resolution of our images is < 1" 
so we can eliminate most background galaxies from the 
GC candidate lists. This careful GC candidate selection, 
and the fact that we can usually trace the spatial distri- 
butions of the GC systems over their full radial extent, 
means that the errors on our derived total numbers and 
specific frequencies are significantly smaller than those 
from past studies, and meet the criteria listed in Sec- 
tion 2.2.1. 

The GC system properties of NGC 1023, NGC 7332, 
NGC 7339, and NGC 7457 were derived from BVR imag- 
ing acquired with the 4096 x 4096-pixel Minimosaic cam- 
era on the WIYN 3.5-m telescope. 2 The GC populations 
of NGC 3379, NGC 3384, NGC 4472, and NGC 5813 
were observed in BVR with the 8192 x 8192-pixel Mo- 
saic imager on the Mayall 4-m telescope at Kitt Peak 
National Observatory. 

Galaxies from the A CS Virgo Cluster Survey. — The GC 
system properties of NGC 4350, NGC 4459, NGC 4473, 
and NGC 4564 come from the ACSVCS survey, which 
was mentioned briefly in Section 2.1. Peng et al. (2008) 
used the ACSVCS data to derive total numbers and spe- 
cific frequencies for 100 early- type galaxies in the Virgo 
Cluster. The ACS FOV is 202" x 202", so the radial 
coverage of the GC system for a galaxy centered in the 
field is 143", which corresponds to 11.5 kpc at 16.5 Mpc, 
their assumed distance to Virgo. In order for a galaxy to 
be included in Table 1, at least 50% of the radial extent 
of the GC system must have been observed. Therefore 
the relevant question is: given the areal coverage of the 
Peng et al. study, for which of the 100 ACSVCS galaxies 
has this criterion been satisfied? 

In Rhode et al. (2010), we used the results from our 
wide-field GC system survey to derive a relationship be- 
tween the stellar mass of a galaxy and the radial extent of 
the galaxy's GC system. Taking that relation (Equation 
1 in Rhode et al. 2010), I calculated the stellar mass of 
a galaxy for which the ACSVCS imaging would encom- 
pass 50% of the radial extent of the GC system. Galaxies 
with log(M/MQ) < 11.3 fall into this category. The gi- 
ant galaxies in the ACSVCS are all elliptical or SO galax- 
ies, so I converted this mass value to an absolute V-band 

2 The WIYN Observatory is a joint facility of the University of 
Wisconsin, Indiana University, Yale University, and the National 
Optical Astronomy Observatory. 
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magnitude by assuming mass-to-light ratios for E and SO 
galaxies from Zepf & Ashman (1993). According to this 
calculation, the ACSVCS imaging should be sufficient to 
meet the 50% radial coverage criteria for ellipticals with 
My > -20.9 and SOs with My > -21.2. 

I searched the list of ACSVCS galaxies to find E/SO 
galaxies that meet these magnitude criteria, have mea- 
sured SMBH masses, and were not already included in 
my sample. The result is a list of six galaxies out of the 
100 ACSVCS targets. Two of these galaxies, NGC 4486 A 
and NGC 4486B, are ellipticals that are companions to 
the giant elliptical NGC 4486 (M87). Peng et al. (2008) 
estimate the total number of GCs in NGC 4486A and 
NGC 4486B is N GC = 11±12 and V GC = 4±11, respec- 
tively. The galaxies therefore fail one of the stated crite- 
ria, that the error on N GC must be < 40%. Four ACSVCS 
galaxies — the two ellipticals NGC 4473 and NGC 4564 
and the two SO galaxies NGC 4459 and NGC 4473 - 
meet all of the requirements and are thus included in 
Table 1. The N GG , My values, and blue GC fractions 
in Table 1 are taken directly from Peng et al. (2008). I 
combined N GG and My with my assumed M/Ly values 
to calculate Sn and T for the table. 

Galaxies from the Gemini Study of Faifer et al. — Two 
of the galaxies listed in Table 1 come from a study by 
Faifer et al. (2011), who used the GMOS camera on 
both the Gemini North and South telescopes to image 
the GC systems of early-type galaxies in the g', r', and 
i' broadband filters. They estimate N GG for two galax- 
ies that have published SMBH masses: NGC 3115 and 
NGC 4649. Based on the relation between galaxy stel- 
lar mass and GC system extent from Rhode et al. (2010), 
the GC system of the SO galaxy NGC 3115 should extend 
to ^26 kpc from the galaxy center, while the elliptical 
galaxy NGC 4649's GC system should go out to a radius 
of ^73 kpc. Faifer et al. trace the radial distribution of 
the GC population of NGC 3115 to -5' (14 kpc) and the 
GCs in NGC 4649 to -8' (40 kpc), so the radial cover- 
age requirement appears to have been met. The galaxies 
were imaged in the g', r', and i' broadband filters and 
Faifer et al. provide values for the fraction of blue GCs 
in the system. 

Faifer et al. fitted both power laws and deVaucouleurs 
profiles to the radial distributions of the GC systems of 
the target galaxies, and integrated the profiles to cal- 
culate Nqc and Sn for two radial limits, 50 kpc and 
100 kpc, presumably because they were uncertain about 
how far out the GC systems extended. To calculate 
Nqc for NGC 3115, I integrated the best-fitting deVau- 
couleurs law profile from Faifer et al. out to 26 kpc (9.3' 
for a distance to the galaxy of 9.7 Mpc). I combined this 
with their assumed absolute magnitude for NGC 3115 
(My = -21.13) and M/L v = 7.6 to calculate S N and T. 
I went through similar steps for NGC 4649, this time inte- 
grating the best-fitting deVaucouleurs GC system profile 
to 73 kpc (15.0' for the assumed distance of 16.8 Mpc) 
and using My = —22.38 and MfLy = 10 to calculate 
Nqc, Sn, and T. The final values for these two galaxies 
appear in Table 1. 

Galaxies from Other Studies in the Literature. — The GC 
system properties of the remaining six galaxies in the 
sample of 20 were drawn directly from various studies in 
the literature that satisfy the selection criteria outlined in 



Section 2.2.1. These galaxies also have published SMBH 
masses. The references from which these quantities are 
taken are listed in the table. Note that the Nqc and 
specific frequency values for the GC system of the Milky 
Way are taken from Ashman & Zepf (1998), who based 
their estimates on mulitiple studies of the Galactic GC 
system. The Nqc and specific frequencies for M31's GC 
population are calculated by combining the results given 
in Ashman & Zepf (1998) and Barmby et al. (2000) (for 
the total number estimates and blue fraction) with Per- 
rett et al. (2002) (for the metal-poor/blue GC fraction). 

2.2.3. Resettling the SMBH Values to Match the 
Distances Assumed for the GC System Studies 

When I assembled the data for Table 1, I noticed that 
the distances to the galaxies that were assumed for the 
GC system studies were often different — sometimes by 
as much as ~10% — from the galaxy distances assumed 
in the studies from which the Msmbh estimates were 
drawn. To make sure that any observed spread in the 
Vcc-AfsMBH points was not due to differences in the 
adopted distance, I have rescaled the SMBH masses to 
the GC system distances (which appear in column 7 of 
Table 1). I followed the general approach of Gultekin et 
al. (2009) and assumed that the Msmbh values are pro- 
portional to the distance, and calculated rescaled Msmbh 
numbers by multiplying them by the ratio of the GC 
distance to the original assumed SMBH distance. The 
SMBH values listed in column 3 of Table 1 are these 
rescaled values. It is worth noting that the scale fac- 
tors applied to the SMBH masses ranged from 0.88—1.08 
(with most in the range 0.95—0.99) and the resultant 
changes to the Msmbh values used here were small and 
well within the stated errors on the Msmbh values given 
in the literature. 

3. RESULTS 
3.1. Relation between Nqq and Msmbh 

The main result discovered by BT10 and investigated 
by HH11 and Snyder et al. (2011) is the correlation be- 
tween the log of Nqc and the log of the SMBH mass for 
early-type galaxies. Figure 1 shows this correlation for 
my sample of giant E, SO, and spiral galaxies. Three of 
the galaxies in the sample — NGC 1399, NGC 3379, and 
NGC 5128 — are plotted twice. This is because they each 
have two SMBH mass determinations that significantly 
differ from one another but are deemed individually re- 
liable in the compilation of Gultekin et al. (2009). Fol- 
lowing that paper and BT10, I show both SMBH masses 
in the plots and have used both estimates in my fits, 
assigning each of the points half- weight. Red filled cir- 
cles denote galaxies with classical bulges and green open 
triangles mark galaxies with bulges that are sometimes 
classified in the literature as possessing a pseudobulge; 
this issue will be discussed in more detail in Section 3.4.2. 

So that I can compare my results directly to those 
of the previous studies on this topic, I have used the 
same statistical analysis methods employed by BT10 and 
HH11. Specifically, I have fitted lines to the data in- 
cluded in Figure 1 using the x 2 -minimization method de- 
scribed in Trcmaine et al. (2002) and Press et al. (1992). 
The method, called f itexy in Press et al. (1992), com- 
putes a best-fit linear relation of the form: 
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Fig. 1. — The log of the number of GCs versus the log of 
the SMBH mass for the 20 giant galaxies included in the 
final sample and listed in Table 1 . The red filled circles 
mark galaxies with classical bulges and the green open 
triangles denote galaxies that may have pseudobulges, 
according to the literature. NGC 1399, NGC 3379, and 
NGC 5128 each have two values of SMBH mass that are 
given half-weight in the fitting process (see Table 1 for 
more information). The solid line is the best- fit line for 
the 18 E/SO galaxies in the sample (Equation 8). The 
dotted line is the best-fit line to the E/SO galaxies from 
BT10 and the dashed line is the best- fit relation for E 
galaxies from HH11. 
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while taking into account the errors in both x and y. The 
quantity 
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is minimized in the fitting process; here, e X i and e V i are 
the errors on x and y, respectively. The errors are as- 
sumed to be symmetric, so I have calculated error bars 
on logarithmic quantities by averaging the high and low 
error bars; i.e., the error on each quantity is: 

criog i = l/2[log(i + CTj) - log(i - <7i)], (5) 

where i in this case is either Nqq or -Msmbh- 

BT10 included thirteen elliptical and SO galaxies in 
their Nqq-M^mbu fit; their best-fit relation is plotted as 
a dotted line in Figure 1 and is given in their paper in 
the form: 
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HH11 fitted the -/Vqc-Msmbh correlation using 18 ellip- 
ticals from their sample of 33 galaxies and termed this 
their "baseline relation" . It is plotted as a dashed line in 
Figure 1 and can be written (following the format of the 
best-fit line from BT10) as: 
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The galaxies with ^(Msmbh/Mq) greater than ^7.5 in 
Figure 1 generally do seem to follow the linear relations 
of BT10 and HH11, but with a few galaxies (NGC 4649, 
NGC 3115, and NGC 4350) deviating from the lines. 
The group of four galaxies with the lowest SMBH masses 
(NGC 7457, NGC7332, the Milky Way, and NGC3384) 
also lie somewhat off the BT10 and HH11 relations. If I 
follow BT10 and include all 18 E/SO galaxies in my fit, 
I derive the following best-fit relation: 
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This relationship is plotted as a solid line in Figure 1. It 
has a steeper slope than the relations from both BT10 
and HH11, although the slope and intercept differ from 
the previous values by less than 3<r, so are nevertheless 
consistent within the errors. 

To assess the quality of the fits in their paper, BT10 
calculated the x 2 per degree of freedom, which is often 
called the "reduced x 2 " ■ Normally when one fits a line, 
the number of degrees of freedom is N ~ 2, where N is 
the number of data points included in the fit. BT10 used 
N — 4 to normalize their x 2 values because three of the 
points they used are related (the three galaxies that each 
have two SMBH mass determinations in Table 1) and the 
inclusion of these points likely reduces the true number of 
degrees of freedom (Tremaine 2011, private communica- 
tion). For their sample of 13 E/SO galaxies, BTlO's best- 
fit A^gc--Msmbh line yields a reduced \ 2 = X 2 /(N - 4) 
= 6.6. For my sample of 18 E/SO galaxies (including 
the same three galaxies with multiple SMBH measure- 
ments) my best-fit relation yields X 1 /{N — 4) = 10.4 and 
X 2 /{N — 2) = 9.3. Therefore, although I have limited my 
sample to galaxies with relatively well-determined iVcc 
values, the result is a Nqc-Msmbh n t 01 lower quality 
(as measured by the reduced x 2 ) than the fit produced 
by BT10 for their galaxy sample. On the other hand, 
the errors on my best-fit slope and intercept in Equa- 
tion 8 are smaller than the errors on those quantities in 
the HH11 fit and comparable to the errors in the BT10 
fit. 

One way that previous authors have used to gauge the 
strength or quality of the correlation, as well as to ac- 
count for the possibility that there is an intrinsic dis- 
persion in the Nqc-Msmbh relation, is to add an addi- 
tional dispersion term to one of the measurement errors 
and repeat the fitting process (cf. Tremaine et al. 2002, 
BT10, Snyder et al. 2011). The assumptions here are 
that the measurement errors are accurate and that there 
is some underlying, real galaxy-to-galaxy scatter present 
in the galaxy properties. With this type of approach, the 
amount of additional scatter that produces a reduced x 2 
~ 1 is inferred to be the true intrinsic dispersion of the 
quantity. To implement this, one replaces the quantity 
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e^j in Equation 4 with (e 2 , + 6 2 ,) 1 / 2 and performs the fit- 
ting procedure, adjusting eo until the value of the reduced 
X 2 is 1. Smaller values of eo are then interpreted as in- 
dicating a tighter correlation between the two quantities 
in the fit. 

For my sample of 18 E/SO galaxies, adding an in- 
trinsic scatter of eo ~ 0.45 dex in quadrature to cither 
log -/Vqc or log Msmbh produces a fit with reduced x 2 ~ 
1. The slope and intercept of the resultant fit are (using 
the notation from Equation 3) /3 = 1.11±0.19 and a = 
8.10±0.11, which are the same within the errors as the 
slope and intercept derived with no additional intrinsic 
scatter (Equation 8). For completeness, I should note 
that another approach sometimes seen in the literature 
is to simply add an intrinsic scatter term directly to the 
measurement errors, i.e., e y i = e y i + eq and then perform 
the fitting process. In this case, eo ~ 0.3 dex (which cor- 
responds to a linear factor of ~2) needs to be added to 
the errors in order to produce reduced x 2 of unity for the 
E/SO galaxy sample. 

When I include all the galaxies in the sample (20 ellip- 
ticals, SOs, and spiral galaxies) and calculate the best-fit 
relation between Agc and Msmbh the result is: 
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Including the two spiral galaxies in the fit steepens the 
slope of the relation not because of the data point for 
M31, which lies very close to the best-fitting relation for 
the E/SO galaxies, but because of the data point for the 
Milky Way, which has comparatively small error bars 
and deviates significantly from the best-fitting E/SO re- 
lation. The quality of this fit (as measured by reduced 
X 2 ) is degraded compared to the fit that includes only 
E/SO galaxies; here, X 2 /(N - 4) = 12.4 and X 2 /(N - 2) 
= 11.2. The scatter as measured by the eo parameter 
is also larger: one must add eo = 0.50 in quadrature to 
the errors in order to produce a reduced x 2 value of ~1. 
Note that the results (slopes, intercepts, x 2 values, etc.) 
of the Nqc— -Msmbh fitting process described in this sec- 
tion and in later sections of the paper are summarized in 
Table 2. The table lists: the sample used in the fitting 
process, the number of points in the sample, the inter- 
cept and slope of the best-fitting relation, the reduced 
X 2 values, and the value of the intrinsic scatter (eo). 

3.2. M-a Relation for the Galaxy Sample 

Another useful or relevant comparison comes from de- 
termining, for my galaxy sample, whether a tighter cor- 
relation exists for Agc-Msmbh or for Msmbh-c- BT10 
found that the former produced a higher-quality fit and 
therefore a stronger correlation. Figure 2 shows the 
Msmbh-c relation for my 20-galaxy sample. As in Fig- 
ure 1, three galaxies are plotted twice and given half- 
weight in the fits. Galaxies with classical bulges are 
again represented by red filled circles and galaxies with 
pseudobulges are plotted with green open triangles; see 
Section 3.4.2 for more discussion. 

I used the f itexy method to fit a line to the 20 galaxies 
in Figure 2. The result is shown as a solid line in the 
figure and can be written: 
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Fig. 2. — The log of the SMBH mass versus the log of the 
velocity dispersion of the galaxy (for early-type galaxies) 
or the bulge (for spiral galaxies) for the 20 giant galaxies 
in Table 1. Red filled circles denote galaxies with clas- 
sical bulges and green open triangles represent galaxies 
classified in the literature as having pseudobulges. Three 
galaxies are plotted twice and given half-weight in the 
fitting process (see Table 1). The solid line is the best- 
fitting line for all 20 galaxies (Equation 10) and the dot- 
ted line is the best-fit M-a relation from BT10. 
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The fit yields reduced x values of x /( N ~ 4 ) = 13.0 
and x 2 / (N — 2) = 11.8. When I limit the galaxy sample 
to only the 18 E/SO galaxies (i.e., excluding the Milky 
Way and M31), the best-fit line becomes: 
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Both the error on the intercept and the reduced x val- 
ues become slightly larger in this case; here, x 2 / (N — 4) 
= 13.8 and \ - 2 /{N - 2) = 12.4. For both the full galaxy 
sample and the E/SO sample, adding an additional scat- 
ter of eo ~ 0.4 dex in quadrature to the errors on the log 
of the SMBH masses brings the reduced x 2 values to ~1. 

The corresponding best-fit M-a relation found by 
BT10 for their sample of 13 E/SO galaxies is plotted as 
a dotted line in Figure 2 and is: 
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BT10 found that both the reduced x 2 value and the 
amount of scatter in their best-fitting Nqq-M^msh re- 
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lation were smaller than those of their best-fitting M-a 
relation; thus, they found that A^gc is a better predictor 
of SMBH mass than is velocity dispersion a. I find mixed 
results when I compare the quality of the A^gc--^SMbh 
and M-a correlations for my galaxy sample. As in BT10, 
the reduced % 2 values of my initial fits are slightly smaller 
for the A"gc--^smbh correlation than for the M-a rela- 
tion. The errors on the coefficients of my best-fit lines 
are likewise smaller for the Agct-^SMbh relation than for 
the M-a relation. However, roughly the same amount of 
scatter (e ~ 0.4 dcx) needs to be added in quadrature 
to the errors on the data points in both the A^gc-^smbh 
and M-a fits in order to bring the reduced \ 2 values to 
unity. 

Finally I should note for the sake of completeness that 
the slope of the SMBH mass — velocity dispersion rela- 
tion derived for my sample of E/SO galaxies (5.34±0.27) 
is consistent within <3 times the errors with the corre- 
sponding relation from BT10, as well as with the derived 
slopes from other recent giant galaxy studies. For exam- 
ple, Ferrarese & Ford (2005) find the slope of the M-a 
relation to be 4.86±0.43, Tremaine et al. (2002) derive 
4.02±0.32, and Gultekin et al. (2009) calculate a slope 
of 3.96±0.42 for ellipticals and 4.24±0.41 for both early- 
and late-types. This is not surprising, of course, since the 
SMBH mass and velocity data I used were largely drawn 
from the same published measurements used by these au- 
thors (although I have rescaled the SMBH masses to the 
distances assumed for the GC system studies). 

3.3. Are Metal-Poor or Metal-Rich GCs Better at 
Tracing SMBH Mass? 

As I mentioned briefly in Section 2.2.1, the broadband 
colors of GCs can be an indicator of their metallicities. 
It is generally true that for old ( > 2—3 Gyr) stellar pop- 
ulations, optical colors primarily trace metallicity, with 
blue colors corresponding to metal-poor populations and 
red to metal-rich. This is caused by a combination of fac- 
tors, such as line blanketing due to the presence of metals 
in stellar atmospheres and metal-rich stars having lower 
equilibrium temperatures on the main sequence and in 
later evolutionary stages (e.g., Schwarzschild et al. 1955, 
Sweigart & Gross 1978). The GC systems of many giant 
galaxies have color distributions with more than one peak 
(e.g., Zepf & Ashman 1993, Kundu & Whitmore 2001, 
Peng et al. 2006), so these multiple peaks seem to imply 
the presence of GC populations with different metallic- 
ities, presumably formed in different episodes. Spectro- 
scopic metallicities, kinematic studies, and near-IR colors 
of GCs confirm this interpretation (e.g., Kundu & Zepf 
2007, Strader et al. 2007, Beasley et al. 2008, Alves-Brito 
et al. 2011). The fact that the Milky Way and M31 also 
have two populations of GCs with different metallicities 
and kinematics (Zinn 1985, Perrett et al. 2002) seems to 
suggest that multiple GC populations may be a common 
characteristic of giant galaxies. 

The overarching goal of many of the studies of ex- 
tragalactic GC systems in recent years has been to un- 
derstand the range of giant galaxy GC system proper- 
ties — including the presence of multiple peaks in GC 
color/metallicity distributions — within the larger con- 
text of galaxy formation and evolution (see the review on 
this subject by Brodie & Strader 2006). Over time a gen- 
eral picture has developed that incorporates aspects of a 



number of different scenarios for GC and galaxy forma- 
tion (e.g., Ashman & Zepf 1992, Forbes et al. 1997, Cote 
et al. 1998, Beasley et al. 2002, Santos 2003, Rhode, Zepf, 
& Santos 2005, Muratov & Gnedin 2010). The precise 
details are not yet well-established, but the data seem to 
point toward a picture (that Brodie & Strader 2006 term 
the "synthesis scenario") in which metal-poor GCs were 
formed when low-mass, gas-rich protogalactic fragments 
collided and merged at high redshift (z^lO— 15) and 
metal-rich GCs were formed during dissipational mergers 
at later times. In this picture the blue, metal-poor GC 
populations in giant galaxies are therefore fossils left over 
from the first epoch of GC and galaxy assembly. The ef- 
fect of cosmological "biasing" may also be an important 
piece of the puzzle; massive galaxy halos located in dense 
environments are expected to have a larger proportion of 
their mass accumulated by a given redshift than lower- 
mass galaxies in lower-density environments (e.g., Shcth 
2003). 

In some models, this initial phase of GC formation 
is truncated or suppressed by some mechanism around 
z^lO for a brief period of time. Santos (2003) suggests 
rcionization would serve as such a truncation mechanism, 
and would suppress both GC and structure formation for 
< 1 Gyr. Forbes et al. (1997) suggest that some more lo- 
calized form of feedback from the first generation of stars 
and GCs could shut off metal-poor GC formation for a 
few Gyr. In the semi-analytic simulations of Muratov & 
Gnedin (2010), blue, metal-poor GCs are formed in early 
mergers of low-mass protogalactic clumps, whereas both 
blue and red GCs are formed in later mergers of more 
massive galaxies. In any case, after the initial epoch of 
galaxy assembly and GC formation, more GCs are ex- 
pected to form each time the host galaxy undergoes a 
gas-rich, dissipational major merger. Because stellar evo- 
lution enriches the interstellar and intcrgalactic medium 
over time, any new GCs that are created in subsequent 
events are expected to be more metal-rich than the ear- 
liest generation. This type of broad outline can explain 
the metallicity gap between GC subpopulations of similar 
age, as well as several other observations, such as: more 
massive galaxies have larger relative numbers of blue, 
metal-poor GCs (Rhode et al. 2005, 2007); the average 
metallicity of metal-poor GCs in more luminous galax- 
ies is slightly higher than that in less luminous galaxies 
(Strader, Brodie, & Forbes 2004); and blue, metal-poor 
GCs have different spatial distributions than red, metal- 
rich GCs in some galaxies (e.g., Zinn 1985, Moore et 
al. 2006, Dirsch et al. 2003, Bassino et al. 2006). On 
the other hand, this rough outline ignores complicat- 
ing effects such as dynamical destruction of GCs over 
time, which may vary depending on the properties of 
the GC itself and/or its parent galaxy (e.g., Vesperini 
2000). It also ignores the possibility that galaxy merg- 
ers may themselves destroy a significant fraction of the 
GCs present in the progenitor merging galaxies, not just 
create new ones (Kruijssen et al. 2012). Both of these 
effects may have substantially modified the trends in GC 
numbers in different galaxies that we measure today (at 
z=0) and may make it considerably more difficult to ac- 
curately interpret our observations. 

A scenario like this may also have implications for the 
observed correlation between A^gc an d Msmbh- BT10 
briefly speculate that the apparent connection between 
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TVqc and -Msmbh may arise because both SMBHs and 
GC populations grow via recent major galaxy mergers. 
Alternatively, SMBHs and GCs may be linked because 
they originate at high-z in gas-rich young galaxies. We 
can consider these two possibilities in the context of the 
broad picture I have outlined of the formation of metal- 
poor and metal-rich GC subpopulations. If SMBHs are 
primordial objects and their properties are set during 
the initial formation epoch of their host galaxies, one 
might expect a strong correlation between the numbers 
of blue, metal-poor GCs in the galaxy and the SMBH 
mass. If, on the other hand, SMBHs are fed primarily 
by recent major galaxy mergers, then one might expect 
a stronger correlation between Msmbh and the numbers 
of red, metal-rich GCs in a galaxy, since the metal-rich 
GC population is likely built up over time during major 
merger events. 

BT10 examined the Nqq— -Msmbh relation separately 
for blue and red GCs in their sample of 13 galaxies but 
found that the correlations were similar in quality; as 
they note, the fractions of blue and red GCs in giant 
galaxies appear to be fairly constant. The blue fraction 
/blue for my sample of 20 galaxies is listed in column 
11 of Table 1. The blue GC fractions for the galaxies 
range from 0.52 to 0.72 and the mean /blue is 0.61, with 
a standard deviation of 0.07. Figure 3 shows the log of 
the number of blue GCs and red GCs versus the log of 
the SMBH mass plotted in the top panel and bottom 
panel, respectively. (Galaxies with classical bulges and 
pseudobulges are again shown with different symbols.) 
The plots look similar and there are only small differences 
in the fitting results from f itexy. Counting only the blue 
GCs, the best-fit relation for the 18 E/S0 galaxies in my 
sample becomes: 
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Fig. 3.— The log of Agc versus the log of Msmbh 
for the blue, metal-poor GCs (top) and red, metal-rich 
GCs (bottom) for the sample of 20 giant galaxies in 
this study. Galaxies with classical bulges and pseudob- 
ulges are marked as in the other figures. NGC 1399, 
NGC 3379, and NGC 5128 each have two values of SMBH 
mass that are given half-weight in the fitting process. 
The solid lines are the best-fit lines derived using the 18 
E/S0 galaxies in the sample (Equations 13 and 14). 
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Including only the red GCs in my E/S0 galaxies yields: 



log 



Ms 



MBH 



Mr, 



.52 ±0.03) + (1.20±0.05)log 



N 



GCrcd 
500 



These relations are shown as solid lines in the appropriate 
plot (blue GCs in the top panel, red GCs in the bottom 
panel) in Figure 3. The blue-only GC sample yields a 
slightly smaller reduced x 2 value than the red GC sam- 
ple; the former yields x 2 /(N-4) = 10.6 for the blue GCs 
vs. x 2 /( N "4) = 11-0 for the red GCs. The blue GC sam- 
ple also needs slightly less additional scatter added to the 
data points than the red GC sample in order to produce 
a reduced x 2 ~ 1; the blue GCs require eo = 0.48 dex vs. 
e = 0.51 dex for the red GCs. Note that the blue GC 
fractions that appear in Table 1 do not have accompany- 
ing uncertainties listed. Assigning an accurate error to 
the blue fraction can be difficult, especially when data are 
being compiled from different studies that use a variety of 
methods for determining blue and red GC fractions (e.g., 
fitting with mixture modeling code, or simply dividing a 
completeness-corrected sample at a certain fiducial color) 
and in many cases no error is given on the estimated blue 
or red fraction. To investigate how possible uncertain- 
ties on the blue and red fraction would contribute to the 



relation between number of blue or red GCs and SMBH 
mass, I incorporated a 10% error on the blue and red 
fractions into the total errors on -/Vcc.biuc and TVccrcd 
and ran f itexy again for the 18 E/S0 galaxies. Includ- 
ing only the blue GCs, and incorporating a 10% error on 
the blue fractions given in Table 1, the best- fit relation 
for the 18 E/S0 galaxies is: 
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Including only the red GCs and including a 10% error on 
the red fraction yields for the E/S0 galaxies: 
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(16) 

With these additional error estimates included, the re- 
duced x 2 f° r the fit to the blue GC sample is x 2 /(N-4) 
= 9.0 and the fit to the red GC sample yields x 2 /(N-4) 
= 9.5. 

If the blue GC population were a better tracer of 
SMBH mass than the red GC population and if the over- 
all picture of hierarchical galaxy and GC formation ac- 
tually holds true, it might suggest that the masses of the 
SMBHs are set during the initial formation epoch of the 
host galaxies and are less dependent on subsequent merg- 
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ers that occur. However, given that reduced x 2 values 
are themselves subject to uncertainty (on the order of a 
few tenths in this case; e.g., Andrae et a!. 2010), the very 
small difference in reduced \ 2 between the blue and red 
GC populations is not significant. 

3.4. Trends with Host Galaxy Morphology 
3.4.1. Splitting the Sample by Hubble Type 

HH11 compiled a larger sample of galaxies than BT10 
with a wider range of morphological types and hence 
were able to investigate how the correlation of A^gc 
with Msmbh varies with galaxy morphology. The HH11 
sample included 21 elliptical galaxies, eight SO galax- 
ies, and four spiral galaxies. My more restricted sam- 
ple of 20 galaxies with reliable GC system properties in- 
cludes 10 ellipticals, eight lenticular galaxies, and two 
spiral galaxies, so I can also explore this issue, although 
to a more limited extent than HH11. Figure 4 shows 
log Nqc versus log Msmbh for the 20 galaxies separated 
by morphological type, with ellipticals, SO galaxies, and 
spiral galaxies plotted from left to right. The best-fit 
ATgc - Msmbh relation for E/SOs (Equation 8) is shown 
as a fiducial; it appears in all three plots as a solid line. 
The three galaxies with two Msmbh values, NGC 1399, 
NGC 3379, and NGC 5128 are (as before) plotted twice; 
all three are Es and thus are included in the leftmost plot. 
The morphological type of NGC 4594 (M104, a.k.a. the 
Sombrero galaxy) is listed in Table 1 as "SO/Sa" and I 
have included this galaxy in the SO plot; it is the point 
with the largest Nqq value, lying very close to the E/S0 
best-fit line. NGC 4594 is often classified by appearance 
as an Sa spiral galaxy, probably at least in part because 
its nearly edge-on orientation makes its dust lane and 
disk appear more prominent. Its quantitative properties, 
however, are more like that of an SO galaxy: its bulge- 
to-disk ratio is ~6, bulge fraction (B/T) is 0.86 (Kent 
1988), and its B - V color is 0.84 (RC3; de Vaucouleurs 
et al. 1991). In any case, it is unclear whether NGC 4594 
belongs in the SO plot (middle) or the spiral galaxy plot 
(right) in Figure 4. 

HH11 found that the elliptical galaxies in their sample 
exhibited a strong correlation between -/Vgc and Msmbh 
and that the spiral galaxies, with the exception of the 
Milky Way, closely followed the same relation. The SO 
galaxies in their sample, on the other hand, showed "no 
clear trend" in a log(JVoc) versus log (M S mbh) plot. The 
results for my sample of 20 galaxies are similar to those 
of HH11. The 10 elliptical galaxies in the left panel of 
Figure 4 closely follow the fiducial E/S0 line. When I 
fit the Ngc— Msmbh relation for just these 10 elliptical 
galaxies, the result is 
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with a reduced \ 2 values of X 2 /(N— 4) = 9.4 and 
X 2 /(N— 2) = 7.7 and eo = 0.37. (For comparison, the 
slope and intercept of the elliptical-only relation from 
HH11 are 0.98±0.10 and 8.30±0.29, respectively; see 
Equation 7.) The slope and intercept of this line are 
the same within the errors as those for my E/S0 sample 
(Equation 8). The spiral galaxy points in the rightmost 
panel in the figure both fall close to the E/S0 line. In the 
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Fig. 4. — The same data as in Figure 1, except with 
the 20 galaxies in the final sample separated by mor- 
phological type. Red filled circles again represent galax- 
ies with classical bulges and green open triangles show 
galaxies with pseudobulges. The solid line is the best-fit 
Nqq — Msmbh relation for E/S0 galaxies (Equation 8). 
NGC 4594 is included in the SO plot, although this galaxy 
is sometimes classified as an Sa. The dotted line in the 
middle panel is the best-fitting relation for the SOs with 
NGC 4594 included, and the dashed line is the best- 
fitting relation with NGC 4594 excluded. 



middle panel of SO galaxies, NGC 4594 (which may be an 
Sa or an SO galaxy) also falls close to the E/S0 best- fit 
line. The other SO galaxies have a relatively small range 
of N GC values (-120-500) but a large range in SMBH 
mass. Fitting a line to the SO points with NGC 4594 
included changes the slope of the best- fit Nqc— Msmbh 
relation to 1.6±0.1. Excluding NGC 4594 and fitting 
the points for the remaining SO galaxies yields only a 
very shallow trend of slightly increasing log(A r Gc) with 
increasing log(AfsMBH) with a slope of 10.2±5.5; this is 
not consistent with the best-fitting relations for the E/S0 
sample or the elliptical-only sample. 

3.4.2. Classical Bulges vs. Pseudobulges 

Before drawing any conclusions about trends in mor- 
phology, it is worth exploring one more morphology- 
related issue and asking whether the behavior of galax- 
ies with classical bulges differs from that of galaxies 
with pseudobulges in the Nqc— Msmbh plane. Classical 
bulges are thought to form via relatively rapid (i.e., evo- 
lution on dynamical time scales), violent processes like 
dissipational collapse and mergers. In contrast, pseudob- 
ulges appear to be built up gradually over a few billion 
years by "secular" , internal processes in galaxy disks, like 
gas transport aided by bar instabilities and spiral struc- 
ture (Kormendy & Kennicutt 2004). By this definition, 
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elliptical galaxies are in the "classical" bulge category 
and SO and spiral galaxies can fall in either category, 
although Kormendy & Kennicutt (2004) note that the 
type of secular evolution that produces pseudobulges is 
most likely to be important in intermediate- to late-type 
galaxies (Sbc galaxies and later). Kormendy, Bender, & 
Cornell (2011) recently analyzed a sample of ^50 spi- 
ral, SO, and elliptical galaxies and found that the disk 
galaxies with pseudobulges do not follow the same M-a 
relation as ellipticals and classical-bulge galaxies. Galax- 
ies in their sample with pseudobulges at least showed a 
much larger scatter in the M-a plane and arguably dis- 
played no correlation at all. 

To examine whether the A^gc —Msmbh correlation 
changes for galaxies with pseudobulges vs. those with 
classical bulges, I searched the literature to find bulge 
classifications for the 20 galaxies in my sample. The 
elliptical galaxies all fall in the "classical" category. I 
found classifications for the SO and spiral galaxies pri- 
marily in the lists included in Kormendy & Kennicutt 
(2004), Fisher & Drory (2008, 2010), and Kormendy et 
al. (2011). Based on my literature search, four of the 20 
galaxies in my sample may have pseudobulges: the Milky 
Way (Shen et al. 2010 and references therein), NGC 3384 
(Kormendy et al. 2011, Pinkney et al. 2003), NGC 7332 
(Pinkney et al. 2003, Falcon-Barroso et al. 2004), and 
NGC 7457 (Kormendy 1993, Kormendy & Kennicutt 
2004). Classifying a galaxy as possessing a pseudobulge 
is notoriously difficult (Kormendy & Kennicutt 2004), so 
none of these four galaxies has been unequivocally shown 
to have a pseudobulge rather than a classical bulge. The 
classifications of NGC 7332 and NGC 7457 are especially 
uncertain. The evidence that NGC 7332 has a pseu- 
dobulge is merely circumstantial: observations with an 
integral-field unit show that this galaxy has a bar and a 
cold, counter-rotating stellar disk, as well as stellar kine- 
matics in its inner regions that are typical of a "boxy 
bulge" (Falcon-Barroso et al. 2004). Together these im- 
ply a pseudobulge rather than a classical bulge. In the 
case of NGC 7457, Kormendy (1993) found that its cen- 
tral velocity dispersion is very small for its bulge luminos- 
ity (i.e., its bulge lies well below the Faber- Jackson rela- 
tion; Faber & Jackson 1976); this is one of the properties 
Kormendy & Kennicutt (2004) identify with pseudob- 
ulges. Pinkney et al. (2003) confirmed' that NGC 7457 
has central stellar kinematics characteristic of a pseudob- 
ulge, but Kormendy et al. (2011) place this galaxy in the 
"classical bulge" category in their paper. 

The four galaxies that may have pseudobulges are des- 
ignated in Figures 1—4 with green open triangles, while 
the classical-bulge and elliptical galaxies are plotted with 
red filled circles. Figure 1 shows the full sample of 20 
galaxies in the Agc~~ Msmbh plane; the four pseudob- 
ulge galaxies are in the lower left of the diagram. In Fig- 
ure 4 the galaxies are split up by morphology, so three 
of the pseudobulge galaxies (NGC 3384, NGC 7332, and 
NGC 7457) appear in the middle panel and the fourth 
(the Milky Way) appears on the spiral galaxy panel on 
the right. Three of the four pseudobulge galaxies are 
systematically above and to the left of the best-fit E/S0 
relation in the Nqq — Msmbh plane, i.e., their GC pop- 
ulations are too large for their SMBH masses (or their 
SMBH masses are too small for their GC numbers). The 
other pseudobulge galaxy (NGC 3384) lies on the rela- 



tion. On the other hand, there are so few pseudobulge 
points and there is so much scatter around the best- 
fitting A^gc— -Msmbh relation that it is not possible to de- 
termine from this sample whether the pseudobulge galax- 
ies are truly outliers or are simply, by coincidence, scat- 
tering in the same direction. 

Lastly, I should note that the galaxies with pseudob- 
ulges seem to fall along the same M-a relation (Fig- 
ure 2) as the rest of the galaxies in the full sample. Two 
of the four galaxies with pseudobulges (NGC 3384 and 
NGC 7332) intersect the best-fitting line for the full sam- 
ple (the solid line in Figure 2), and a third galaxy (the 
Milky Way) lies very close to the line. Only NGC 7457 
scatters far away from the best-fit relation (and its classi- 
fication as a pseudobulge galaxy was uncertain) . As men- 
tioned earlier, Kormendy et al. (2011) concluded that 
pseudobulge galaxies do not follow the M-a relation for 
ellipticals and classical-bulge galaxies. The galaxies that 
seemed to deviate from their best-fitting M-a relation in 
Kormendy et al. were located in a specific region of the 
M-a plane: below SMBH masses log(M/M ) <8 and 
velocity dispersions log(cr) <2.2. Kormendy et al. (2011) 
had ~15 galaxies with pseudobulges in that area of the 
plane. My sample has very few objects in that region of 
the M-a plane so it is difficult to assess whether the data 
presented here contradict their conclusions or not. 

3.4.3. What Might the Trends in Morphology Mean? 

In part because of the small numbers of SO, spiral, 
and pseudobulge galaxies in the sample, it is not clear 
what the pseudobulge vs. classical bulge values in the 
Nqq— Msmbh plane and the data in Figure 4 are indi- 
cating. It is possible that there is: (1) a marked trend 
of increasing Aqc with increasing Msmbh mass for el- 
liptical galaxies and perhaps spiral galaxies, but a much 
weaker trend for SOs; or (2) a clear trend in increasing 
Nqc with Msmbh for classical-bulge galaxies but not for 
pseudobulge galaxies; or (3) an overall trend in Nqq vs. 
Msmbh for giant galaxies, but with plenty of scatter, 
especially at lower A^qc and lower SMBH masses. 

The situation implied by the first possibility does not 
have an obvious physical explanation; it is not clear what 
mechanism would cause SO galaxies to have a limited 
range of Aqc values but a wide range of SMBH masses 
(log(MsMBH/M Q ) <~ 6—9) while ellipticals have Aqc and 
Msmbh values that are closely linked. Whatever direct 
or indirect correlation might be present between GC pop- 
ulations and SMBH masses in elliptical galaxies may, for 
some unknown reason, break down for most SO galaxies. 

The second possibility, on the other hand, could make 
sense in terms of how we think both SMBHs and bulges 
grow. Kormendy et al. (2011) succinctly describe two 
different "modes" for the feeding and growth of SMBHs, 
that echo the two different growth modes for building 
galaxy bulges. One mode is violent, rapid mergers that 
drive gas infall and efficiently build up both classical 
bulges and the seed SMBHs within them. The other 
mode is the more gradual, local, stochastic growth (e.g., 
build-up of a central concentration in a disk galaxy via 
disk instabilities and gas-flow along bars) that builds 
pseudobulges, and that may also result in slower and 
more modest growth of SMBHs. If GC populations, clas- 
sical bulges, and SMBHs can all be grown and signif- 
icantly increased via mergers, one might imagine that 
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galaxies with classical bulges are then also likely to have 
Nqc values and Msmbh values that increase together. 
Galaxies with more quiescent histories, that perhaps 
made GCs in their early assembly stages at high-z but 
did not do much more to supplement their GC popu- 
lations over time, perhaps would end up today as disk 
galaxies with a pseudobulge, a modest-sized SMBH, and 
a weaker correlation between their Msmbh and Ncc val- 
ues. In other words, whatever correlation between Nqq 
and Msmbh was set up in the earliest stages of a galaxy's 
assembly may be strengthened in galaxies that continue 
to grow via major mergers, but weakened in galaxies with 
pscudobulges that evolve more quiescently. 

Finally, the third scenario mentioned above is that the 
Nqq~ Msmbh correlation is generally present in giant 
galaxies but is weaker and has more scatter at lower 
Nqq values and lower SMBH masses, and by implica- 
tion, lower-mass galaxies, since both Ncc and Msmbh 
correlate with galaxy mass (see next section) . If this is 
the case, the explanation could be similar to the sce- 
nario for possibility (2). Perhaps a common mechanism 
like major galaxy mergers drives the formation of very 
massive SMBHs and very populous GC systems in mas- 
sive galaxies and this mechanism becomes less important 
in the evolution of lower-mass galaxies, their GC popu- 
lations, and their SMBHs. The net effect would be that 
the GC numbers and SMBHs "decouple" (because the 
phenomenon that linked them in the first place, major 
gas-rich mergers, is no longer of primary importance in 
their evolution) and the correlation gets weaker with de- 
creasing galaxy mass. 

3.5. Correlations of Nqc with Other Galaxy Properties 

3.5.1. Number of GCs Versus Galaxy Stellar 
Luminosity and Mass 

The number of GCs in some giant galaxies has been 
observed to follow a rough correlation with galaxy lu- 
minosity, with Ncc K L a , where a is between 1 and 2 
(Ashman & Zepf 1998 and references therein). A sen- 
sible next step is to use the current galaxy sample to 
examine this and other correlations that Nqq has with 
galaxy properties, as well as to compare the strength 
of this correlation to the Ncc~ Msmbh relation. To do 
this, I have supplemented the galaxies in the N=20 sam- 
ple that I used in earlier sections of the paper with ten 
more galaxies with measured Ncc from my wide-field 
GC system survey These ten additional galaxies have 
wcll-dctcrmincd GC system properties (i.e., the criteria 
in Section 2.2.1 are met), but the galaxies were not in- 
cluded in the previous analysis because they do not have 
measured SMBH masses. The properties of these addi- 
tional ten galaxies are listed in Table 3. The columns 
in the table are: galaxy name and morphological type; 
predicted SMBH mass calculated from the best-fitting 
Nqc— Msmbh relation (Equation 9; see Section 3.6 for 
details); central velocity dispersion; My for the galaxy; 
distance to the galaxy in Mpc; number of GCs; GC spe- 
cific frequencies Sn and T; fraction of blue GCs in the 
GC system; and the reference for the GC system prop- 
erties. 

The left-hand side of Figure 5 shows the log of Ncc 
versus the log of the y-band luminosity in solar units, 
Ly/L Ql for the 30 galaxies listed in Tables 1 and 3. The 



right-hand side of the same figure shows log of Agc plot- 
ted against the log of the galaxy stellar mass in solar 
masses (M*/M ) for the same 30 galaxies. To calculate 
stellar masses for the galaxies, I combined their My val- 
ues with the Zepf & Ashman (1993) mass-to-light ratios 
discussed in Section 2.2.1. In both plots, galaxies from 
Table 1 with classical bulges are designated with red filled 
circles and galaxies with pscudobulges with green open 
triangles, while the 10 additional galaxies from the wide- 
field survey are plotted with open stars. I searched the 
literature for bulge classifications for these additional 10 
galaxies and found that one of the objects, the Sb spi- 
ral galaxy NGC 891, may have a pseudobulge; de Jong 
et al. (2008) characterize its bulge as a boxy pseudob- 
ulge based on surface photometry. The point represent- 
ing NGC 891 is labeled in the plots. The point with 
the largest Ncc error bars in both plots is NGC 7331; 
this galaxy has a relatively low inclination (i~75 degrees) 
compared to the other spiral galaxies in the survey, which 
makes it difficult to detect GCs against the background 
of the galaxy's spiral disk. The end result was a Ncc 
value that was relatively uncertain (Rhode et al. 2007). 
Lastly, note that there are only 30 points in the figure 
because the three galaxies that were plotted twice in the 
figures that included SMBH mass are only plotted once 
here (they each have only one value of Ncc and galaxy 
luminosity or stellar mass). 

Fitting the 30 data points shown in the panels of Fig- 
ure 5 yields the following best-fit linear relations: 

Tat - Nc^r -1 
log-^ = (10.44 ± 0.01) + (0.70 ± 0.02)log-^ (18) 

and 

lo s£r = ( 1L33 ± °- 01 ) + (°- 78 ± °- 02 ) lo g^T ( 19 ) 

M 500 

Written another way, these relationships mean Ncc °^ 
(Ly/L©) 1 - 42 and Ncc oc (M^/Mq) 1 - 28 , respectively. 
The reduced \ 2 value for the fit between Ncc and to- 
tal galaxy stellar mass is xV(-^ — 2) = 8.9, compared 
to 15.8 for the Ncc~ luminosity fit. This makes sense 
given the appearance of the data: the log of the number 
of GCs shows a tight correlation with the log of the host 
galaxy stellar mass, whereas the log(Acc) vs. log(Ly) 
data show much more scatter, especially for the lower- 
luminosity galaxies with log(Ly) < 10.75. The intrinsic 
scatter is substantially smaller in this relation compared 
to any of the Nqq— Msmbh relations: adding e = 0.14 
dex in quadrature to the errors brings the reduced % 2 
value to ~1. (For completeness, I also checked the fit- 
ting result for the 20-galaxy sample in Table 1, with- 
out the additional ten galaxies from the survey. For the 
N=20 sample, the intrinsic scatter remains small, at eo 
= 0.15 dex.) The correlation between Ncc and galaxy 
stellar mass in Figure 5 is the strongest of any of the cor- 
relations analyzed in this paper. It is important to note 
that the reduced % 2 and eo values would be even smaller 
if errors on the galaxy luminosity or mass were included 
in the fits. Specifically arbitrarily setting the galaxy 
mass error to 10% of the stellar mass yields a best-fitting 
line with the same slope and intercept (within the un- 
certainties on those coefficients) as the fit without errors 
on the mass, but with a x 2 /(A — 2) = 5.1. Setting the 
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Fig. 5. — The log of Nqc versus the log of the l/-band 
luminosity and galaxy stellar mass for the 20 galaxies in 
Table 1 and the ten galaxies from my wide-field survey in 
Table 3. Galaxies with classical bulges and pseudobulges 
in the Table 1 sample are marked with red filled circles 
and green open triangles, respectively, as in the other fig- 
ures. Galaxies from the wide- field survey are designated 
with open stars. The spiral galaxy NGC 891 (labeled on 
the plots) is not shown in green but may have a pseu- 
dobulge. The best-fit lines are plotted on the figures and 
given in Equations 18 and 19. The relation between Nqc 
and total galaxy luminosity yields \ 2 1 (N — 2) = 15.8 and 
the relation between -/Vqc and total galaxy stellar mass 
yields X 2 /{N - 2) = 8.9. 



relative error on the galaxy stellar mass to 25% yields the 
same coefficients and y^/{N — 2) = 1.5. Even without 
error bars in the ^-direction, the quality of the correla- 
tion between Nqq and galaxy mass is much better than 



the Ngc — Msmbh and M-a correlations explored in the 
earlier sections of the paper. 

When I include only the blue, metal-poor or red, metal- 
rich GCs in Nqc and fit the data points versus total 
galaxy stellar mass, I find that the slope of the line re- 
mains the same within the errors, but that the quality of 
the fit is better for the blue GCs than for the red ones. 
The former yields % 2 /(iV - 2) = 8.7 and e = 0.14 dex, 
whereas the latter yields ^/(N — 2) = 11.1 and eo = 
0.15 dex. The slightly improved fit yielded by the blue 
GCs may suggest that it is chiefly the blue, metal-poor 
GCs that are driving the correlation between Nqq and 
galaxy stellar mass. 

One object in Figure 5, the SO galaxy NGC 7457, has 
a particularly large value of Nqq for its stellar mass, and 
lies significantly off the best-fit relation. It is the lowest- 
mass galaxy in the sample, with log(M*/M Q ) = 10.6 but 
it has an estimated population of 210±30 GCs (Hargis 
et al. 2011). This galaxy is in a low-density environment 
(p = 0.13 Mpc" 3 ; Tully 1988) and its GC system was 
studied as part of our wide-field survey. There is not an 
obvious reason for its high Nqc value and high specific 
frequency (S N = 3.1±0.7). Excluding NGC 7457 from 
the sample and fitting log(iVGc) vs. log(M*/M Q ) yields 
the same coefficients as in Equation 19 within the errors 
and improves the fit, so that x 2 /{N — 2) becomes 6.0 and 
eo becomes 0.10 dex. 

Lastly, I should note that in addition to calculating to- 
tal galaxy stellar masses from the y-band luminosity, I 
also computed stellar masses from if -band luminosity of 
each galaxy and used those to examine the relationship 
between Nqc and galaxy stellar mass. I combined the 
total if-band apparent magnitude listed in the NASA 
Extragalactic Database (NED) — which in turn come 
from either from the 2MASS Extended Source Catalog 
(Jarrett et al. 2000) or the 2MASS Large Galaxy Atlas 
(Jarrett et al. 2003) — with the adopted distance modu- 
lus to derive total if -band luminosity for each galaxy. I 
then applied the appropriate if -band mass-to-light ratio 
from Spitler et al. (2008) for each galaxy's morphological 
type to calculate galaxy stellar mass. I fitted log(iVGc) 
vs. log(M*/M Q ) for these if -band-derived masses and 
found that the best-fitting coefficients were similar to 
those in Equation 19 but the scatter in the plot was larger 
and the reduced \ 2 values increased by a factor of ~1.5. 
So at least for this specific data set, the F-band galaxy 
magnitudes and stellar masses produced better-quality 
results in terms of the correlation between number of 
GCs and galaxy stellar mass. 

3.5.2. Number of GCs Versus Bulge Luminosity and 

Mass 

The analysis in the previous section examined corre- 
lations of Nqq with the total stellar mass and luminos- 
ity of the host galaxy, but there are reasons why look- 
ing at Nqq versus the mass and luminosity of just the 
spheroidal component might make more sense. As men- 
tioned in Section 1, the analysis of Snyder et al. (2011) 
indicated that the correlation between Nqc and AisMBH 
arises because both of those quantities correlate with 
the binding energy or potential well depth of the bulge 
component of the host galaxy. It has been argued that 
at least some GCs are actually associated with galaxy 
bulges rather than with thick disks or halos. Specifically, 
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the properties of the more metal-rich GC subpopulation 
of the Milky Way may indicate it is a bulge population 
(e.g., van den Bergh 1993). In this type of picture, the 
red, metal-rich GCs would then be an analogous popu- 
lation in other spiral and elliptical galaxies (e.g., Forbes, 
Brodie & Larsen 2001). In previous extragalactic GC 
system studies, the specific frequency Sn is sometimes 
normalized using bulge luminosity rather than total lu- 
minosity for spiral galaxies, (e.g., Harris 1981, Cote et al. 
2000, Forbes et al. 2001). For these reasons, I have also 
investigated the relationship between JVqc and galaxy lu- 
minosity and mass using the bulge luminosity and mass 
of the spiral galaxies rather than the total stellar lumi- 
nosity and mass. 

To calculate bulge luminosities and masses for the sam- 
ple galaxies, I found published estimates of the mean 
bulge fraction (B/T, the fraction of the total luminos- 
ity that comes from the bulge) for galaxies of various 
Hubble types. I used the data given in Binney & Merri- 
field (1998), which is based on the work of Kent (1985), 
and combined it with the total l/-band magnitudes and 
galaxy stellar masses in Tables 1 and 3 to calculate My 
and log(M/M ) for each galaxy's bulge. I assumed bulge 
fractions of 1.0 for the elliptical galaxies in the sam- 
ple, 0.65 for the lenticular galaxies, 0.375 for the Sab 
galaxies, 0.25 for the Sb galaxies, 0.2 for the Sbc galax- 
ies, and 0.15 for the Sc galaxy. For NGC 4594, I used 
the B/T measurement from Kent (1988) that was men- 
tioned in Section 3.4.1. I also tried using bulge frac- 
tions derived from data given in Simien & de Vaucouleurs 
(1986), which vary from the Binney & Merrifield values 
by a few tenths, and finally I tried applying the bulge-to- 
disk ratios for spiral galaxies given in Graham & Worley 
(2008). For both the Simien & de Vaucouleurs (1986) 
bulge-to-disk ratios and the Graham & Worley (2008) 
ratios, the resultant -/Vqc vs. bulge luminosity and bulge 
mass fits yielded similar slopes and intercepts but pro- 
duced slightly more scatter and higher % 2 values. I have 
adopted the Binney & Merrifield bulge fractions as the 
final values. 

The left panel in Figure 6 shows the log of Nqq versus 
the log of the y-band luminosity of the host galaxy bulge 
for the 30 galaxies in Tables 1 and 3. The symbols are 
the same as in Figure 5. The positions of the data points 
for the elliptical galaxies have not changed compared to 
Figure 5, but the spiral and SO galaxies (in the lower- 
luminosity region of the diagram) show smaller scatter 
and a more linear behavior. For example, NGC 7457 
was the galaxy with the lowest luminosity in Figure 5 
and was an outlier, whereas in Figure 6 it lies within the 
cluster of other data points at the low-luminosity end of 
the sample. The right-hand panel in Figure 6 shows log 
of iVcc plotted against the mass of the galaxy bulge com- 
ponent. Some of the outliers in Figure 5 — specifically, 
NGC 7457 and NGC 7331 — also have become less dis- 
crepant in this figure compared to the plot of log(iVGc) 
versus total galaxy mass. 

The quantitative results from the fitting process are 
mixed. The log of Nqc correlates much more closely with 
the log of the bulge luminosity than it does with total 
galaxy luminosity. However, the correlation of log(7Vcc) 
with bulge mass is much weaker than that with total 
galaxy stellar mass. The best-fitting lines to the 30 data 
points in Figure 6 are 
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Fig. 6. — The log of Nqc versus the log of the 1^-band 
luminosity or stellar mass of the bulge for spiral galax- 
ies or the entire galaxy for elliptical galaxies. The data 
points represent the 30 galaxies in Tables 1 and 3. The 
symbols are the same as in Figure 5. The spiral galaxy 
NGC 891 is not shown in green but may have a pseudob- 
ulge. The best-fitting relations are shown as a dashed 
lines and given in Equations 20 and 21. The relation 
between Nqc and bulge luminosity gives x 2 /(N — 2) = 
9.2 and the relation between Nqq and bulge mass gives 
X 2 /(N-2) = 11.5. 
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log -= = (11.11 ± 0.01) + (1.15 ± 0.03)log^ 

(21) 

which translate to the proportionalities Nqq oc 
(iv^buigc/i©) 1 - 02 and A GC cx (M* Mlge / M Q ) - 87 . The 
fit of log(iVGc) vs. bulge luminosity yields "x 2 /{N — 2) = 
9.2 and e = 0.17 dex, whereas the linear fit of log(AGc) 
vs. bulge mass gives x 1 / (N — 2) = 11.5 and e = 0.23 dex. 
By these measures, the tightest correlation of log(iVGc) 
with any of the galaxy properties examined here is still 
the correlation with galaxy stellar mass, shown in the 
right-hand panel of Figure 5 and given in Equation 19. 

One might assume that if it is the red, metal-rich GCs 
that are associated with a galaxy's spheroidal compo- 
nent, then perhaps the numbers of red GCs will better 
correlate with bulge luminosity or mass. Including only 
the red GCs in log(iVGc) an d fitting the log of the bulge 
luminosity or mass produces linear relations with the 
same slope (within the errors) as those in Equations 20 
and 21, but larger \ 2 and eo values. In fact, I also tried 
executing the fits counting only the blue GCs in each 
galaxy and found that counting the blue GCs produces a 
better-quality fit than counting only the red ones. Still, 
using the total number of GCs produces the best results 
in terms of x 2 and eo- 

To summarize, the relationship between -/Vgc and 
galaxy mass, or even with bulge luminosity or mass, 
for the 30 galaxies in this sample yields reduced x 2 val- 
ues that are smaller than or comparable to the small- 
est reduced x 2 values yielded by the -/Vgc —Msmbh fits 
in Section 3.1. Furthermore, the intrinsic scatter is 
much smaller for the fits explored in this section (eo ~ 
0.1—0.2 dex) than for either the Nqq— Msmbh or the 
M— a relations (both of which have eo -0.4 dex). This 
seems to imply that the relationship between Nqc and 
the luminosity or mass of the host galaxy is the more 
fundamental one, and that the observed correlation be- 
tween -/Vgc and SMBH mass results from both of these 
quantities being correlated with galaxy mass. As noted 
earlier, Snyder et al. (2011) came to a similar conclusion 
based on their analysis of the residuals in the relation be- 
tween -/Vgc and Msmbh for a sample of elliptical galaxies 
drawn from BT10 and HH11. The result here is there- 
fore not unexpected since — as the previous papers on 
this topic (BT10, HH11, and Snyder et al. 2011) have 
all noted — it seems highly unlikely that the GCs in a 
galaxy will have some sort of direct causal link to the 
SMBH located in the galaxy center. 

It is relevant to discuss here a recent paper by Sadoun 
& Colin (2012) that examined other correlations between 
GC systems and SMBH masses. For a sample of twelve 
giant spiral, SO, and elliptical galaxies, Sadoun & Colin 
(2012) find a correlation between the projected velocity 
dispersion of the GC system (ogc) and the SBMH mass 
that is as tight as the M-a relation for those twelve galax- 
ies. In the discussion section of their paper, Sadoun & 
Colin point out that their observational results might be 
understood in the context of recent numerical simulations 
of galaxy formation done by Jahnke & Maccio (2011). 
Jahnke & Maccio simulate the hierarchical assembly and 
merger history of dark matter halos over cosmic time, 
with added recipes for star formation, black hole accre- 



tion, and bulge evolution in the simulated galaxies. They 
then examine the scaling relation between Msmbh and 
bulge mass that is produced in the resultant galaxies 
at z—0. Their conclusion is that correlations between 
SMBH mass and bulge mass (or velocity dispersion as 
a proxy for mass) do not require a direct physical link, 
but rather can arise through normal galaxy merging in 
a ACDM universe that occurs over the course of a giant 
galaxy's history, from high-redshift until z—0. In Jahnke 
& Maccio 's simulations, the scaling relations are a nat- 
ural outcome of the simultaneous build-up of both the 
central black hole and the galaxy's stellar component by 
major galaxy mergers. Jahnke & Maccio state that, "we 
in principle expect a correlation with Msmbh for any 
(mass) parameter that is subject to the same ACDM as- 
sembly chain", including correlations between GC pop- 
ulations and Msmbh- Sadoun & Colin (2012) point out 
that this could provide an explanation for their observed 
Msmbh and ctgc results, and likewise it seems to provide 
a broad context within which we might understand the 
correlations between the various parameters presented in 
the current paper. 

3.6. Deriving an M-a Relation for the Wide-Field GC 
System Survey Galaxy Sample 

One last piece of analysis that can be done with the 
data set compiled here is to calculate predicted SMBH 
masses for the giant spiral and elliptical galaxies in Ta- 
ble 3. These ten galaxies have well- measured Nqc val- 
ues from my wide-field survey but do not yet have mea- 
sured masses in the literature for the SMBH that pre- 
sumably exists in each of their central regions. Nine of 
the ten galaxies in Table 3 are spiral galaxies and one 
(NGC 4406) is a massive Virgo Cluster elliptical. Conse- 
quently it seems appropriate to derive Msmbh from the 
relationship between Nqq and Msmbh given in Equa- 
tion 9, i.e., the linear relation based on the full sample of 
elliptical, SO, and spiral galaxies analyzed in Section 3.1 
and listed in Table 1, rather than using the relation de- 
rived from only the E/S0 galaxies (Equation 8). The cal- 
culated SMBH masses are listed in column (2) of Table 3. 
The errors on the masses are computed by determining 
the predicted SMBH mass for the upper and lower end 
of the possible Nqq values. For example, Nqc for the 
spiral galaxy NGC 1055 is 220±30, so the uncertainties 
in Msmbh reflect the possible range of SMBH values if 
A GC is 220-30=190 or 220+30=250. Published mea- 
surements of the central velocity dispersion for nine of 
the ten galaxies are given in the third column of the ta- 
ble. One galaxy (NGC 7339) apparently does not have a 
published velocity dispersion, so its SMBH mass is listed 
in the table but it is excluded from the rest of the analysis 
in this section. 

Figure 7 shows the nine galaxies from Table 3 with pre- 
dicted SMBH masses and measured velocity dispersions, 
added to the M — a plot shown earlier in this paper (Fig- 
ure 2). As before, the red filled circles are classical-bulge 
galaxies and the green open triangles are pseudobulge 
galaxies. The open stars are the nine new galaxies from 
the wide-field GC system survey. The two lines are the 
same as those shown in the first M — a plot: the solid line 
is the best-fitting M — a relation given in Equation 10 
and the dotted line is the best-fitting M — a relation 
published in BT10. 
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Fig. 7. — The log of the measured or predicted SMBH 
mass versus the log of the central velocity dispersion for 
the 30 giant galaxies in Tables 1 and 3. As in previ- 
ous figures, galaxies from Table 1 are plotted with red 
filled circles (classical bulges) or green open triangles 
(pseudobulges) and the Msmbh data plotted are mea- 
sured values. Open stars denote galaxies from Table 3 
and the Msmbh values are predicted based on the num- 
ber of GCs in the galaxy and Equation 9. The lines are 
the same as those plotted in the earlier version of this 
figure, Figure 2: the solid line is the best-fitting line for 
the 20 galaxies in Table 1 and the dotted line is the M-a 
relation from BT10. 



Most of the galaxies from the wide-field survey do seem 
to follow the general relationship between Msmbh and 
velocity dispersion shown by the two lines. If one adds 
the data points for these nine galaxies to the data points 
from Table 1 and fits a linear relation to the 29 points, 
the result is a line with the same slope and intercept 
(5.88±0.19 and 8.19±0.03, respectively) within the er- 
rors as those of Equation 10, and similar reduced x 2 val- 
ues (x7(A - 4) = 12.5 and x*/(N - 2) = 11.6). On 
the other hand, a few spiral galaxies from the survey — 
namely, NGC 1055, NGC 3556, and to a lesser extent 
NGC 4013 — deviate noticeably from the relations in 
that they have very large Msmbh values for their mea- 
sured velocity dispersions. One relevant question to ask 
is: do these galaxies have large predicted SMBH masses 
because they have anomalously large GC populations? 
The answer is no: all three of these galaxies have the ex- 
pected number of GCs based on their stellar mass, i.e., 
they are not outliers on the Agc versus mass plot (see 
Figure 5 and the data in Table 3). What the data seem 
to be indicating instead is that Agc is not a particularly 
good predictor of SMBH mass for some galaxies. The 
three galaxies that deviate most strongly in the M — a 
plot are all spiral galaxies. HH11 pointed out that the to- 
tal number of GCs in the Milky Way was also not a good 



predictor of its SMBH mass; that is, it deviated strongly 
from their best-fitting Agc~ Msmbh relation and Fig- 
ure 1 shows that it deviates from the relation derived 
here. 

The one additional elliptical galaxy from the wide-field 
survey is NGC 4406, which appears at Zo^A/bh/Mq) ~ 
9.1 and log(cr) ~ 2.4 in Figure 7. The spiral galaxies 
all have relatively low Msmbh and velocity dispersion 
values and lie in the lower left region of the figure. If 
one considers only the spiral galaxies from the survey, it 
appears that they actually do not seem to follow an M — a 
relation and it is really the data point for NGC 4406 
that is driving the agreement between the new data and 
the original best-fitting M — a line. In fact, excluding 
NGC 4406 from the sample and fitting a line to the other 
eight data points from the wide-field survey yields a line 
with zero slope within the errors: 



log- 
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(7.19 ± 0.23) + (0.17 ± 0.80)log 7 



'200 km s- 1 
(22) 

Although it is not entirely obvious what this means, what 
it seems to suggest is that, again, Agc nray not be an 
accurate predictor of SMBH mass, at least not for spiral 
galaxies; if there were a tight correlation between the two 
quantities, one might expect the spiral galaxies to follow 
the expected M—a relation more closely. A larger sample 
of galaxies with well-determined values of Agc, Msmbh, 
bulge velocity dispersion, and bulge classification would 
presumably help to clarify these issues. 

4. SUMMARY 

I have assembled a sample of 20 giant galaxies with 
both measured SMBH masses and reliable estimates of 
Aqc and GC color fractions based on high-quality wide- 
field imaging data. The sample includes eight galaxies 
from my ongoing GC system survey, four galaxies from 
the ACSVCS survey (Peng et al. 2008), two galaxies from 
a Gemini study by Faifer et al. (2011), and six galaxies 
from various studies in the literature. Half of the galaxies 
in the sample are ellipticals, eight are SO galaxies, and 
two are spiral galaxies (the Milky Way and M31). Four 
of the galaxies may possess pseudobulges according to 
previous studies of their light profiles and/or kinematics. 
The sample of 20 galaxies is used to explore correlations 
between the galaxies' GC populations and the masses of 
their central SMBHs. The main findings of this investi- 
gation are as follows: 

1. The E/S0 galaxies in the sample follow a relation 
between the log of Agc and the log of the SMBH 
mass of the form: 



log 



Msmbh 
M m 



.04 ±0.03) + (1.22 ±0.06)log 



N G c 
500 



(23) 

When the two spiral galaxies are included in the 
sample, the slope and intercept of the relation be- 
come 1.52±0.06 and 7.91±0.03, respectively. The 
coefficients of the best-fitting relation for E/S0 
galaxies agree within the errors with those of the 
corresponding relation for early-type galaxies from 
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BT10 (who found a slope of 1.08±0.04 and an in- 
tercept of 8.14±0.04) and the relation for ellipti- 
cal galaxies from HH11 (who derived a slope of 
0.98±0.10 and an intercept of 8.30±0.29). 

2. The M-a relation for the 20 galaxies in the sam- 
ple agrees with those in the recent literature within 
the errors. Three of the four galaxies with pseudob- 
ulges (including the Milky Way) closely follow the 
best-fitting linear relation defined by the classical- 
bulge galaxies. 

3. The numbers of blue, metal-poor GCs in the galax- 
ies yield a slightly smaller reduced \ 2 than the num- 
bers of red, metal-rich GCs, although this small dif- 
ference is probably not statistically significant. In 
the current picture of galaxy and GC system for- 
mation, metal-poor GCs originate in the earliest 
epoch of galaxy assembly. If metal-poor GC popu- 
lations were actually more tightly correlated to the 
SMBH masses in the host galaxies, this might im- 
ply that the correlation between Nqq and Msmbh 
is put in place early in the history of the galaxy 
and does not depend strongly on the occurrence of 
major merger events later. 

4. When the galaxy sample is divided according to 
Hubble type, the elliptical galaxies retain the tight 
correlation between A^gc and SMBH mass seen in 
the full sample, whereas the SO galaxies show more 
scatter. This is similar to the result in HH11, who 
found that the SO galaxies show a large disper- 
sion in black hole mass for a given Nqq- The four 
galaxies that may have pseudobulges also seem to 
show increased scatter compared to the elliptical 
and classical-bulge galaxies. However, the number 
of pscudobulge galaxies in the sample is small so it 
is not clear whether or not the latter trend is real 
and therefore implies something about how SMBHs 
grow. 

5. Ten more galaxies from my wide-field GC system 
survey were used to supplement the original sample 
and explore correlations between Nqq and other 
galaxy properties. The strongest correlation with 
the smallest intrinsic scatter is a correlation be- 
tween the number of GCs and the total stellar 



mass of the host galaxy. In general the correla- 
tions between Nqq an d galaxy mass and luminos- 
ity or bulge mass and luminosity are much tighter 
than any of the Nqq — -Msmbh correlations for this 
sample. This seems to confirm the idea that the 
observed connection between A^gc and Msmbh in 
giant galaxies is a consequence of the connection 
between both of these quantities and the galaxy 
potential. 



6. Finally, the Agc - Msmbh relation derived here 
is used to calculate predicted SMBH masses for 
the ten additional galaxies with measured Agc 
from the wide-field survey but without existing 
Msmbh measurements in the literature. The single 
elliptical galaxy in this subsample lies close to the 
expected M — a relation, but the spiral galaxies 
show larger scatter, suggesting that Agc is n °t a 
reliable predictor of SMBH mass for some galaxies. 
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TABLE 1 

Giant Galaxies with Measured SMBH Masses and Global GC System Properties 



Name Type M BH Ref a M v Dist N G c S N T f hluc Ref 

(M ) (km/s) (Mpc) 



N821 


E6 


3.7 


+2.5, 


-0.7 


xlO 7 


12 


209±10 a 


-21 





22.4 


320+45 


1 


3+0 


2 


1 


5+0 


2 


0.70 


22 


N1399t 


Elpec 


1.1 

4.5 


+0.4, 
+0.6, 


-0.6 
-0.6 


xlO 9 
xlO 8 


12 


337±16 a 


-22 


3 


18.5 


5800+700 6 


7 


0+0 


8 


8 


2+1 





0.60 


3,5 


N3379t 


El 


1.1 
3.6 


+0.7, 
+0.9, 


-0.5 
-0.9 


xlO 8 
xlO 8 


7,12 


206±10 a 


-20 


9 


10.6 


274+31 


1 


2+0 


3 


1 


4+0 


3 


0.70 


20 


N4374 


El 


1.6 


+ 1.2, 


-0.6 


xlO 9 


12 


296±14 a 


-22 


1 


18.4 


1775+381 


2 


7+0 


6 


3 


1+0 


7 


0.59 


10 


N4472 


E2 


1.8 


+0.6, 


-0.6 


xlO 9 


8 


310±10 b 


-23 


1 


16.7 


5870+680 


3 


6+0 


6 


4 


2+0 


6 


0.60 


19 


N4473 


E5 


1.2 


+0.4, 


-0.8 


xlO 8 


12 


190±9 a 


-20 


7 


15.3 


376+97 


2 


0+0 


5 


2 


3+0 


6 


0.57 


17 


N4564 


E6 


6.4 


+0.4, 


-0.9 


xlO 7 


12 


162±8 a 


-19 


9 


15.8 


213+31 


2 


2+0 


3 


2 


6+0 


4 


0.62 


17 


N4649 


E2 


4.8 


+ 1.1, 


-1.1 


xlO 9 


21 


385±19 a 


-22 


4 


16.8 


4073+820 


4 


5+0 


9 


5 


2+1 


1 


0.52 


6 


N5128+ 


Epcc 


3.1 
7.3 


+0.4, 
+ 1.4, 


-0.2 
-3.9 


xlO 8 
xlO 7 


12 


150±7 a 


-22 





4.6 


980+120 


1 


5+0 


2 


1 


8+0 


2 


0.53 


15 


N5813 


El-2 


7.0 


+ 1.1, 


-1.1 


xlO 8 


11 


237±3 b 


-22 


3 


32.2 


3100+800 


3 


9+0 


9 


4 


5+1 


1 


0.57 


13 


N1023 


SBO 


4.2 


+0.5, 


-0.5 


xlO 7 


11 


205±4 b 


-21 


1 


11.4 


494+80 


1 


7+0 


3 


2 


7+0 


4 


0.58 


23 


N3115 


SO 


9.1 


+5.1, 


-2.8 


xlO 8 


12 


230±ll a 


-21 


1 


9.7 


439+130 


1 


6+0 


3 


2 


4+0 


7 


0.52 


6 


N3384 


SBO 


1.8 


+0.1, 


-0.3 


xlO 7 


12 


143±7 a 


-20 


5 


11.6 


122+30 





8+0 


2 


1 


2+0 


4 


0.60 


13 


N4350 


SO 


5.5 


+4.0, 


-4.0 


xlO 8 


11 


181±7 b 


-20 


1 


16.5 


196+60 


1 


8+0 


5 


2 


7+0 


8 


0.71 


17 


N4459 


SO 


7.0 


+ 1.3, 


-1.3 


xlO 7 


12 


167±8 a 


-20 


9 


16.1 


218+28 


1 


0+0 


1 


1 


5+0 


2 


0.52 


17 


N4594 


SO/Sa 


6.6 


+0.4, 


-0.4 


xlO 8 


16 


238±41 c 


-22 


4 


9.8 


1890+200 


2 


1+0 


3 


3 


2+0 


5 


0.63 


20 


N7332 


SOpec 


1.3 


+0.6, 


-0.5 


xlO 7 


12 


125±3 b 


-20 


8 


23.0 


175+56 





9+0 


3 


1 


4+0 


4 


0.57 


23 


N7457 


SO 


3.9 


+ 1.1, 


-1.6 


xlO 6 


12 


69±4 b 


-19 


5 


13.2 


210+30 


3 


1+0 


7 


4 


8+1 


1 


0.58 


14 


Milky Way 


Sbc 


4.3 


+0.4, 


-0.4 


xlO 6 


9 


116±2 d 


-21 


3 


0.0083 


160+20 





5+0 


1 


1 


1+0 


1 


0.70 


1 


M31 


Sb 


1.4 


+0.9, 


-0.3 


xlO 8 


4 


170±5 b 


-21 


8 


0.8 


450+100 





9+0 


2 


1 


6+0 


4 


0.72 


1,2,18 



References. — (1) Ashman & Zepf 1998; (2) Barmby ot al. 2000; (3) Bassino ct al. (2006); (4) Bonder ot al. 2005; (5) Dirsch ot al. 
2003; (6) Faifer et al. 2011; (7) Gebhardt et al. 2000a; (8) Gebhardt ct al. 2012, in preparation; (9) Gillesson ct al. 2009; (10) Gomez & 
Richtler 2004; (11) Graham 2008; (12) Gultekin et al. 2009; (13) Hargis & Rhode 2012; (14) Hargis et al. 2011; (15) Harris et al. 2004; 
(16) Jardel ct al. 2011; (17) Peng ct al. 2008; (18) Pcrrctt et al. (2002); (19) Rhode & Zepf 2001; (20) Rhode & Zepf 2004; (21) Shen & 
Gebhardt 2010; (22) Spitler et al. 2008; (23) Young, Dowcll, & Rhode 2012. 

NOTE. — Morphological types arc from the NASA/IPAC Extragalactic Database (NED), with three exceptions. NED lists NGC 5128 
as an "SO pec" or "E/SO", but studies of this galaxy's properties and stellar populations (e.g., Hui ct al. 1993, Peng ct al. 2004, Harris 
et al. 2004) describe it as an elliptical. NGC 4594 is often classified as an Sa but it has a bulge fraction of B/T— 0.86 (Kent 1988) and 
broadband colors like an SO, so I list it here as "SO/Sa". The Milky Way is listed as an Sbc, which is the classification commonly used 
in the literature (e.g., Ashman & Zepf 1998). 

^Two estimates of SMBH mass arc given for the galaxies NGC 1399 and NGC 5128 in Gultekin ct al. (2009), as well as for NGC 3379 
in BT10. I adopt the approach taken by Gultekin ct al. (2009) and BT10 and use both SMBH mass measurements for these galaxies, 
but increase the error bars by a factor of v^2, to give each measurement half-weight. 

a Vclocity dispersion from Gultekin et al. (2009). 

b Vclocity dispersion from HYPERLEDA. 

c Vclocity dispersion from Jardel ct al. (2011). 

^ Velocity dispersion from Howard ct al. (2008). 

e The GC system properties of NGC 1399 were calculated by Spitler et al. (2008), using measurements drawn from both Dirsch ct al. 
(2003) and Bassino ct al. (2006). 



TABLE 2 

Results of Fitting Process for Nqc and -Msmbh Data 



Sample 


-^points 


a 


P 


X 2 /(N - 4) 


X 2 /(N~2) 




E/SOs 


21 


8.04+0.03 


1.22+0.06 


10.4 


9.3 


0.45 


All 


23 


7.91+0.03 


1.52+0.06 


12.4 


11.2 


0.50 


E/SOs, blue GCs 


21 


8.32+0.03 


1.23+0.06 


10.6 


9.5 


0.48 


E/SOs, red GCs 


21 


8.52+0.03 


1.20+0.05 


11.0 


9.8 


0.51 


Es 


13 


8.17+0.04 


0.95+0.07 


9.4 


7.7 


0.37 


SOs (with N4594) 


8 


7.98+0.05 


1.63+0.14 


16.2 


10.8 


0.70 


SOs (without N4594) 


7 


10.62+2.52 


10.21+5.46 


9.8 


5.9 


0.75 



Note. — For the fits between TVqc and Msmbh, the form of the linear relation is log — ^ BH — a 



+ /31og 5 *?p . Values of eo in column (7) are derived by adding eo m quadrature to the errors on the 
data points. 



Globular Cluster-Supermassive Black Hole Relation 



TABLE 3 

Giant Galaxies with Measured GC System Properties from the Wide-Field Survey and Predicted 

SMBH Masses 



Name 


Type 


Predicted M BH 


a 


M v 


Distance 


^GC 




T 


/blue 


Ref 








(km/s) 




(Mpc) 












N4406 


E3 


1.2(+0.2,-0.2)xl0 9 


235±3 b 


-22.3 


16.7 


2900±400 


3.5±0.5 


4.1±0.6 


0.62 


3 


N891 


Sb 


4.1(+1.9,-1.6)xl0 6 


73±10 b 


-20.8 


8.4 


70±20 


0.3±0.1 


0.6±0.3 


0.50 


1 


N1055 


Sb 


2.3(+0.5,-0.5)xl0 7 


80±15 b 


-21.0 


16.3 


220±30 


0.9±0.2 


1.7±0.5 


0.59 


5 


N2683 


Sb 


9.2(+5.1,-4.2)xl0 6 


118±9 b 


-20.5 


7.7 


120±40 


0.8±0.4 


1.4±0.7 


0.68 


4 


N3556 


Sc 


3.5(+1.6,-1.4)xl0 7 


79±10 a 


-21.2 


12.4 


290±80 


0.9±0.4 


2.2±0.9 


0.64 


4 


N4013 


Sb 


1.2(+0.3,-0.2)xl0 7 


86±9 a 


-20.4 


15.1 


140±20 


1.0±0.2 


1.9±0.5 


0.74 


1 


N4157 


Sb 


5.0(+2.0,-1.8)xl0 6 


90±4 b 


-20.4 


14.7 


80±20 


0.6±0.3 


1.1±0.6 


0.55 


4 


N7331 


Sb 


2.2(+2.3,-1.7)xl0 7 


137±4 b 


-21.7 


13.1 


210±130 


0.5±0.3 


0.8±0.7 


0.50 


4 


N7339 


Sbc 


4.5(+0.9,-0.9)xl0 6 




-20.4 


22.4 


75±10 


0.5±0.2 


1.5±0.5 


0.60 


5 


N7814 


Sab 


1.5(+0.4,-0.3)xl0 7 


170±8 b 


-20.4 


13.2 


165±25 


1.3±0.4 


2.2±0.8 


0.58 


2 



References. — (1) Rhode, Windschitl, & Young 2010; (2) Rhode & Zepf 2003; (3) Rhode & Zopf 2004; (4) Rhode et 
al. 2007; (5) Young, Dowell, & Rhode 2012. 



Note. — Morphological types are from the NASA/IPAC Extragalactic Database (NED). 

a Velocity dispersion from Ho et al. (2009). 
b Velocity dispersion from HYPERLEDA. 



